NASA CFM52126 


'"tf»ASA-CR-T52126) A MACH LINE . PANEL. METHOD ’ ,N 7 8- 27087^ 

FOR COMPUTING THE' LINEARIZED -SUPERSONIC, FLOW 
OVER PLANAR WINGS (Boeiiig Commercial 

Airplane Co. e Seattle) 91 p HC.A0.5/MF A-01 Unclas 

, . ’ CS.CL :0TA G-3/0’2 ' 25873 


A Mach Line Panel Method for 
Computing the Linearized Supersonic Flow 
Over Planar Wings 


F. Edward Ehiers and Paul E. Rubbert 

Boeing Commercial Airplane Company 
Seattle, Washington 


Prepared for 

NASA Ames Research Center 
under Contract NAS 2-7729 



National Aeronautics and 
Space Administration 



1978 




. TECHNICAL REPORT STANDARD TITLE PAGE 


1. Report No 

NASA CR-152126 

2. Government Accession No. 

3. Recipient's Catalog No. 

4 Title and Subtitle 


5 Report Date 

A Mach Line Panel Method for Computing the 

May 1978 

Linearized Supersonic Flow Over Planar Wings 

6. Performing Organization Code 

7 Author(s) 


8. Performing Organization Report No. 

F. E. Ehlers and P. E. Rubbert 


D6-46373 


9. Performing Organization Name and Address 

The Boeing Commercial Airplane Company 

P. O. Box 3707 

Seattle, Washington 98124 


1 2 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 


15 Supplementary Notes 


16. Abstract 

A method is described for solving the linearized supersonic flow over planar wings using panels 
bounded by two families of Mach lines. Polynomial distributions of source and doublet strength 
lead to simple, closed-form solutions for the aerodynamic influence coefficients, and a nearly 
triangular matrix yields rapid solutions for the’ singularity parameters. 

The source method was found to be accurate and stable both for analysis and design boundary 
conditions. Similar results were obtained with the doublet method for analysis boundary 
conditions on the portion of the wing downstream of the supersonic leading edge, but 
instabilities in the solution occurred for the region containing a portion of the subsonic leading 
edge. Research on the method was discontinued before this difficulty was resolved. 


1 7. Key Words 

Linearized theory 
Subsonic flow 
Panel methods 
Aerodynamics 


18. Distribution Statement 

Unclassified-Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No of Pages 

22. Price 

Unclassified 

Unclassified 


86 



10. Work Unit No 


11. Contract or Grant No 

NAS2-7729 

13. Type of FI eport and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


Form DOT F 1700.7 (8-691 



















FIGURES 


No. Page 

1 Mach Line Coordinates x,y 63 

2 Mach Line Paneling on Planar Wing at M =y/2~ 64 

3 Analytic Continuation of the Doublet Distribution on the Wing Associated 

With an Interior Mach Line Panel at M = \fT~ 65 

4 Range of Interior Mach Line Panel and Analytic Continuation of Doublet - 

Strength Outside of the Panel 66 

5 Domain of Dependence of an Interior Panel With Illustration of Upstream 

Characteristic Strips 67 

6 Control Point Locations and Continuity Conditions for Doublet Mach Line 

Paneling 68 

7 Analytic Continuation of the Doublet Distribution on the Wing Associated 

With a Mach Line Panel Adjacent to Supersonic Leading Edge at M =y/2 .... 69 

8 Range of Integration for Supersonic Leading Edge Mach Line Panel Analytic 

Continuation of Doublet Strength Outside of the Panel 70 

9 Range of Integration for Subsonic Leading Edge Mach Line Panel and Analytic 

Continuation of Doublet Strength Outside of the Panel 71 

10 Pressure Distributions on Portion of the Wing Behind Supersonic Leading Edge . . 72 

1 1 Pressure Distributions on Portion of the Wing Behind Subsonic Leading Edge 

Along the Line y = 1.0603 . 73 

12 Error in Downwash From Mach Line Doublet Panel Method Behind 

Subsonic Leading Edge Along the Line y = 1.0603 74 

13 Discontinuity in Pressure Coefficient Along Panel Edges of Mach Line 

Doublet Panels in the Region Behind Subsonic Leading Edge . 75 

14 Pressure Distribution Along Line L Downstream of Special Mach Line With 

Flat Plate Solution Superimposed on Doublet Panel Solution 76 

1 5 Pressure Distribution Along Line L Downstream of Special Mach Line With 

Flat Plate Solution Superimposed on Doublet Panel Solution 77 

1 6 Pressure Distribution Along Line L Downstream of Special Mach Line With 

Flat Plate Solution Superimposed on Doublet Panel Solution 78 

17 Pressure Distribution Along Line L Downstream of Special Mach Line With' 

Flat Plate Solution Superimposed on Doublet Panel Solution 79 

18 Control Point Locations and Continuity Conditions for Source Mach Line 

Paneling 80 

1 9 Comparison of Pressure Distribution From the Source Mach Line Panel Method 

With the Exact Solution From Linearized Theory for Wing With Parabolic Arc 
Profile 81 

20 Comparison of Source Distribution With Wing Slopes 82 

21 Comparison of Downwash From Supersonic Source Design Panel Method 

With Actual Wing Slopes ' 83 

22 Schematic Paneling on Wing Used to Test Mach Line Panel Methods 84 

23 Region of Integration for Leading Edge Source Panel With the Zone of 

Dependence Terminated by a Subsonic Leading Edge 85 

24 Region of Integration for Interior Source Panel With the Zone of Dependence 

Terminated by a Subsonic Leading Edge 86 


I 


v 



CONTENTS 


Page 

1.0 SUMMARY * 1 

2.0 INTRODUCTION 2 

3.0 SYMBOLS AND ABBREVIATIONS 3 

4.0 PLANAR MACH LINE DOUBLET PANEL METHOD WITH 

LINEARIZED BOUNDARY CONDITIONS 5 

4.1 Theoretical Description of the Method 5 

4.2 Discussion of Results From the Planar Mach Line Doublet Panel Method . . 10 

5.0 PLANAR MACH LINE SOURCE PANEL METHOD WITH 

LINEARIZED BOUNDARY CONDITIONS 12 

5.1 Theoretical Description of the Method 12 

5.2 Discussion of Results From the Planar Mach Line Source Panel Method . . . 13 

6.0 CONCLUSION ' 15 

APPENDIX A— DERIVATION OF BASIC FORMULAS FOR THE PLANAR MACH 1 6 

LINE PANEL METHODS. 

A1 Formulas for the Downwash from Supersonic Leading Edge 

Doublet Panels 16 

A2 Formulas for the Downwash from Interior Doublet Panels 18 

A3 Formulas for the Downwash from Subsonic Leading Edge 

Doublet Panels 19 

A4 Description of Panel and Parameter Numbering System 20 

A5 The Doublet Distribution for Panels in Columns Behind the 

Supersonic Leading Edge (i < i c ) 20 

A6 Doublet Distribution for Panels in Columns Behind the Subsonic 

Leading Edge(i >-i c ) 23 

A7 Formulas for the Downwash in Columns Behind Supersonic 

Leading Edge(i < ij 26 

A8 Formulas for the Downwash in Columns Behind Subsonic 

Leading Edgefi > i c ) 27 

A9 Matrix Equations for the Doublet Panel Method 29 

A10 Source Distribution for Columns Behind Supersonic Leading Edge 32 

A1 1 Pressure Coefficient from the Supersonic Leading Edge Panels and . 

from the Interior Panels 35 

A12 Pressure Coefficient from the Subsonic Leading Edge Panels 39 

A 13 Correction to the Pressure Coefficient from Supersonic Leading Edge 
Panels and from Interior Panels When the Subsonic Leading Edge 
Panel of the Same Row is in the Zone of Dependence of the Control Point . 40 
' A14 Source Distribution for Panels in Columns Behind the Subsonic 

Leading Edge 42 


iii 



CONTENTS (Concluded) 

Page 

A1 5 Complete Pressure* Co efficient from the Source Panels 45 

APPENDIX B— INTEGRATION OF THE AERODYNAMIC INFLUENCE COEFFICIENTS 

FOR THE PLANAR MACH LINE PANELS 48 

B 1 Basic Integrals for Supersonic Leading Edge Panels 48 

B2 Basic Integrals for Subsonic Leading Edge Panels 49 

B3 Aerodynamic Influence Coefficients for Supersonic Leading Edge 

Doublet Panels 50 

B4 Aerodynamic Influence Coefficients for Interior Doublet Panels 53 

B5 Aerodynamic Influence Coefficients for Subsonic Leading Edge 

Doublet Panels 54 

B6 Aerodynamic Influence Coefficients for Supersonic Leading Edge 

Source Panels 56 

B7 Aerodynamic Influence Coefficients for Rectangular Interior 

Source Panels 58 

B8 Aerodynamic Influence Coefficients for Subsonic Leading Edge 

Source Panels 59 

B9 Subsonic Leading Edge Correction to the Aerodynamic Influence 
Coefficients for Supersonic Leading Edge Source Panels and 
Interior Source Panels 61 


iv 



A MACH LINE PANEL METHOD FOR COMPUTING THE 
LINEARIZED SUPERSONIC FLOW OVER PLANAR WINGS 


F.E. Ehlers and P.E. Rubbert 
Boeing Commercial Airplane Company 


1.0 SUMMARY 


This report describes a study that was conducted in an attempt to develop advanced super- 
sonic panel methods in linearized supersonic flow. This study centered around the use of 
panels bounded by the two families of Mach lines, a feature which appeared at the outset to 
offer significant advantages. Mach line paneling allows the accurate treatment of discontinui- 
ties in velocity and velocity gradient which occur along special Mach lines emanating from 
discontinuities in the wing planform. The use of Mach line coordinates also leads to simple, 
closed-form solutions for the aerodynamic influence coefficients. By ordering the panels 
along characteristic or Mach line strips a nearly triangular matrix results for the simultaneous 
equations to be solved for the flow parameters. The solution of the flow over the wing is, 
therefore, very rapid. 

Both doublet and source panel formulations were implemented and evaluated. The doublet 
method was tested on a swept, parabolic-cambered wing, and the source method was tested 
on a swept wing with a symmetrical parabolic arc profile. The source method was accurate 
and stable both for analysis boundary conditions in which wing slope is prescribed and for 
design boundary conditions in which surface pressure is prescribed and wing slope is computed. 
Refining of the panel size improved the accuracy of the source method. Similar results for 
analysis boundary conditions with the doublet method were obtained for the portion of the 
wing bounded on the downstream side by two Mach lines which intersect at an interior point 
and inclose only the supersonic leading edge. However, instabilities in the solution occurred 
for the region containing a portion of the subsonic leading edge. 

Research on the method was discontinued before the difficulty was resolved in favor of a 
method (ref. 1) that appeared to offer more promise in terms of its adaptability to general 
nonplanar configurations. This latter method subsequently underwent extensive develop- 
ment and will be reported in a forthcoming NASA contract report. 



2.0 INTRODUCTION 


Since wings in supersonic flow are usually thin and adequately described by linearized boundary 
conditions, it is natural to represent the wing by doublet and source distributions on a plane. 

In contract report NASA CR-2423, Mercer, Weber, and Lesferd (ref. 2) proposed a method 
utilizing Mach line paneling on a plane with linearized boundary conditions and presented 
some preliminary analysis and classification of panel forms which occur in a representative 
wing planform. This report describes the continued investigation into this method. 

The approach appears to have several advantages: 

1. The Mach lines are well defined on the mean plane on which thickness and camber can 
be described by source and doublet distributions. 

2. As discontinuities in velocity and velocity gradient occur only across certain Mach 
lines, denoted here as “special Mach lines,” the discontinuities can be conveniently 
taken into account in the influence coefficients by appropriate paneling. 

3. With the two families of Mach lines as coordinates, the region of integration for the 
source and doublet integrals on the plane is easily described and leads to simple, closed- 
form expressions for the influence coefficients. 

4. Solutions can be obtained extremely rapidly since the aerodynamic influence matrix 
can be placed in nearly triangular form by ordering the panels along characteristic 
strips on the wing planform. 


In this report both the source and the doublet Mach line panel methods are described in 
detail. The source method for the thickness problem is very efficient and yields accurate 
results for both analysis and design type boundary conditions. The doublet method for 
the camber problem was successful in regions directly behind the supersonic leading edges 
but exhibited instabilities in regions downstream of subsonic leading edges. However, because 
^oftKelsuccess of the subsonic panel method of Johnson and Rubbert (ref. 1) and its greater 
applicability to a wide range of configurations and its apparent' extendability to supersonic 
flow, the Mach line panel method described here was abandoned before the difficulty with 
the doublet method in the region behind subsonic leading edges could be resolved. 


The method discussed here is based on linearized supersonic flow theory. The perturbation 
velocity potential 4 is given by the Prandtl Glauert differential equation 


(l - M 2 V 


X X ^v V z ^ 

A O A 0 y o y o ^0*0 


where M is the freestream Mach number and Xq, yQ, zq are the Cartesian coordinates in 
the flow. The freestream is in the Xq direction and the wing boundary conditions are 
applied on the zq = 0 plane. 
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3.0 SYMBOLS AND ABBREVIATIONS 


a,b,c,d Coefficients in polynomial distribution of doublet strength 



C(x,y) Correction function for Mach line strips with a portion of the subsonic . 

leading edge in zone of dependence (sec. B9 of app. B) 

Cjj, Contribution to pressure coefficient from the kth panel in the ith column 

Cp Pressure coefficient 

fj(x>y) Basic functions for integrals over leading edge panel i = 1 to 3 for supersonic 

leading edge panels and i = 4 to 6 for subsonic leading edge panels 

G(i,k,x,y) Polynomial used in doublet strength for the kth panel in the ith column 

gj(x,y) Downwash functions associated with the polynomial coefficients of doublet 

. strength from supersonic leading edge panels; i = 1 ,2,3,4 

gj(x,y) Pressure functions associated with source strength polynomial coefficients 

for supersonic leading edge; j = 1,2, 3, 4 

hj(x,y) Downwash functions associated with the polynomial coefficients of doublet 

strength from interior panels; i = 1,2, 3, 4 

h(x,y) Pressure function associated with the single coefficient of interior source panel 

Hq = \/(x-|)(y-n) Hyperbolic distance on the plane z = 0 in Mach line coordinates 

iorj Column index 

i c Column index for last column behind supersonic leading edge 

k Index for panel number in any column 

M Freestream Mach number 

N(i,k j) Doublet parameter numbers for the kth panel in the ith column; j = 1 ,2,3,4 

N p (i,k) = i(i-l)/2 +k, for i < i c ; panel number for the kth panel in the ith column; 

Np(i,k) = [i c (i c +l) + (3i c - i + 2 ) (i - i c -1)] /2+k, i > i c 

Sj Width of column in region behind supersonic leading edge 

S 2 Width of column in region behind subsonic leading edge 
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SjCx.y) 

Sj(x,y) 

x 0 > y 0’ z 0 

x n 

x,y,z 

x,y 

x i 

y j 

y j 

w 

e l 

e 2 

M 

^ik 

a 

°ik 

4> 


Downwash functions associated with the coefficients of the doublet strength 
in subsonic leading edge panels; i = 1,2, 3 ,4 

Pressure function associated with the coefficients -of the source strength for 
subsonic leading edge panels; j = 1,2 

Cartesian coordinates;-,Xg aligned with the freestream 

Doublet parameters consisting of the polynomial coefficients of the double 
strength for all panels 

Mach line coordinate system, defined -in equation (3) 

Denote differentiation with respect to the designated variable when used as 
a subscript 

x + i e s j ; panel coordinate for panels in the ith row from y coordinate axis 
= y - jsf, j <i c ; panel-coordinate in jth" column for j < i c 
— y - i c Sj - (j - i c ’+ 1) S 2 ? panel coordinate for jth column for j > i c 
Downwash 

Slope of supersonic leading edge defined by x + ejy = 0 
Slope of subsonic leading edge defined by x - e 2 Y - 0 
Doublet strength 

Total doublet strength in kth panel of ith column; equations (A27) and (A34) 
Source strength 

Total source strength in kth panel of the ith column; equation (A64) 
Perturbation velocity potential 
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4.0 PLANAR MACH LINE DOUBLET PANEL METHOD WITH 
LINEARIZED BOUNDARY CONDITIONS 


4. 1 THEORETICAL DESCRIPTION OF THE METHOD 


The camber of a three dimensional wing of small curvature can be simulated by a doublet 
distribution on a plane. In supersonic linearized theory, the perturbation velocity potential 
0 induced by such a doublet distribution is given by the integral 


<K x o> Vo* z o) 2 tt 9z 0 fJ 

s w 


P(£>h) d£dp 



0 ) 


where S w is the po rtion o f the wing cut by the upstream Mach cone from the point 
x 0 ,y 0 ,z (J and B =-\/M2-l. The velocity potential is normalized to the freestream velocity 
and the coordinates x,y,z to the wing chord. The doublet distribution n is determined 
by satisfying the downwash boundaiy conditions on- the plane Zq — 0 

a0/.az o = f Xo (x o ,y o ) ■ (2) 


where f Xq (x 0 ,y 0 ) is the slope of the camber surface zq = f (xq,yq) . The Mach line panel 
method of solving the boundary value problem of equation (2) is facilitated by introducing 
Mach line coordinates defined by the transformation (see fig. 1 ). 


z 0 = z 

x 0 - ByQ = 2Bx/M (3) 

xq + ByQ = 2By/M 


Since the Jacobian of the transformation in equation (3) is 2B/M^, the potential in equation 
(1) in the new variables takes the form 


* 2?rM dz ff 
S w 

Under the transformation in equation (3), the differential equation for 0 becomes 

^zz “ ^ 0 X y = 0 

Differentiating the potential with respect to z and using the fact that the integral as well 
as 0 itself satisfies the differential equation, we can replace the second derivative with 
respect to z by the second derivative with respect to x and y. With this change in 
derivatives, the evaluation of 0 Z on the plane z = 0 may be found by settings = 0 in the 
integrand without first performing the differentiation (ref. 3). Thus we obtain 
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( 4 ) 


_M _j£_ Y r X) d£dp 

2ir 3x9y J J V (x-z)(y-7?) 

where | = £(jj) and y = y(x) are the equations of the leading edge. 

In this form the integration must be performed before differentiation because of the inverse 
square root singularity of the integrand. This problem is eliminated by introducing the new 
variables of integration 

£'~\/x-£ V '~Vy-v ( 5 

which lead to 

\/y- y(x) y/x - i(y - 7? 2 ) 

J (6: 

where we have dropped the primes for convenience. Some simplification results from per- 
forming the differentiation before integration. When equation (6) is substituted into the 
boundary condition in equation (2), with the quantity fx 0 ( x o=y 0 ) accordingly transformed 
into the Mach line coordinates in equation (3), an integral equation for the doublet distribu- 
tion ju(x,y) results. 

Once the doublet distribution on the plane has been solved, the pressure on the wing can be 
evaluated. From the relationship implicit in equation (1), 

fiim 0 = - m/2 
z-K) + 

the pressure coefficient on the upper side of the wing is 


M 

C p =~2 0 Xq = ~’g" + 

(7) 

= (tb) + ^y) 

The pressure coefficient on the lower side of the wing has the same magnitude, but opposite 
in sign. 


The numerical approach to solving the integral equation using the Mach line panel method 
is to divide the wing planform into small panels by a grid system of the two families of Mach 
lines with x and y being constant, and then determine the doublet distribution on each panel 
by imposing the appropriate boundary conditions at certain discrete points on the planform. 
A typical paneled wing with M = yfq. is shown in figure 2. The paneling is comprised of 
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parallelogram-shaped interior panels having edges bounded by two families of Mach lines 
and triangular panels adjacent to the leading and trailing edges. The doublet distribution 
on each panel is approximated by a polynomial with sufficient boundary and continuity 
conditions to establish the coefficients of the polynomial. The doublet strength must vanish 
along a leading and side edge and be continuous across panel edges. The downwash at 
certain points on the panels computed by summing the integrals of the form of equation (6) 
over all panels in the region of dependence must equal to the prescribed camber slope at 
these same points. 

It is necessary that the doublet distribution be continuous over the planform as such dis- 
continuities at panel edges introduce free vortex lines. For convenience of the discussion, 
we shall draw the paneling and wing planform in the oblique Mach line coordinate x,y system 
as if it were orthogonal (which is the case for a Mach number of 1 / 2 ). The interior panels 
are then rectangular. 

In supersonic flow a rectangular panel influences the flow field in the region within the 
downstream Mach cone which emanates from the upstream vertex of the panel. With Mach 
line paneling, the panel edges coincide with the boundaries of this region as seen in figure 3. 
Continuity of the doublet strength is assured by choosing a polynomial which vanishes on 
the lines x = constant and y = constant defining the two upstream boundaries of the panel. 

In terms of local panel coordinates with the origin at the upstream vertex of the rectangular 
panel, the doublet strength p must be proportional to x,y for rectangular panels within the 
planform; i.e., 


H~xy P(x,y) 

where P(x,y) is a polynomial. If we further postulate that on the downstream panel edges 
that the doublet strength be quadratic, this restricts g to the form 

p = xy (a + bx + cy + dxy) ' (8) 

The relation for ju is continued analytically along characteristic strips downstream of the 
panel trailing edges by its quadratic relation on the downstream boundaries: x = Ax and y = Ay; 
where Ax and Ay are the panel length and width, respectively. Hence the contribution to ju in the 
two downstream characteristic strips containing the panel is quadratic in the Mach line variable 
running normal to the edge (see fig. 3). In the remaining portion of the domain of influence 
bounded by the two characteristics or Mach lines from the downstream vertex, the doublet 
strength’ ju is constant as indicated in figure 3. The range of integration for control points 
in the four distinct regions of influence of the panel is illustrated in figure-4 with the corres- 
ponding functional relations for the doublet strength described in equation (8). 

The complete doublet distribution on a given interior panel consists of the function in 
equation (8) associated with that panel plus contributions from: 

1. The quadratic doublet strengths at the downstream edges of those panels in the two 
characteristic strips running upstream of the panel as shown in figure 5. 
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2. The constant doublet strengths at the downstream comers of those panels within the 
region between the two upstream running characteristic strips. 

Equation (8) for the doublet strength has four free parameters, which can be determined 
as follows. 

Continuity of the normal derivative (in which Mach line variables) is applied at the down- 
stream ends of the leading edges of the panel, leaving two parameters to be determined by 
the downwash boundary conditions at'two control points on the panel, see figure 6. How- 
ever, for an interior panel in which one of the leading edges is a special Mach line along which 
discontinuities in velocity or velocity gradient may occur, the requirement of the continuity 
of the doublet gradient is relaxed on this edge and is replaced by a downwash boundary 
condition. The formulas for the downwash induced by a doublet distribution on interior 
panels are derived in section A2 of appendix A. They are simple, rational algebraic expres- 
sions. 

The treatment for a triangular panel adjacent to a supersonic leading edge is similar. Since 
the doublet strength must vanish on the leading edge of the wing, the doublet distribution 
must be of the form 


p = (x + e^y) (a + bx + cy + dxy) (9) 

where x + e^y - 0 defines the portion of the leading edge covered by the panel and the 
origin is on the upstream corner of the panel illustrated in figure 7. It is easily seen that 
p vanishes on the x and y Mach line boundaries of the zone of influence for the leading 
edge panel. In the same manner as for the interior panels, the value of doublet strength is 
continued along the two characteristic strips by the quadratic relations on its downstream 
boundaries. In the region defined by the two x and y Mach lines from the downstream 
vertex the doublet strength p is constant. The range of integration for control points in 
the four distinct regions of influence of the panel is illustrated in figure 8. The four free ' 
parameters in equation (9) are determined as follows. 

In order to ensure the continuity of the doublet gradient along all straight supersonic leading 
edges, we require the polynomial (a + bx + cy + dxy) at the upstream corner to agree with 
the adjacent leading edge panel value as indicated in figure 6. The downwash boundary 
conditions at the upstream corner point, the middle of the panel leading edge, and near the 
downstream corner point furnish the-rest of the conditions to render the four parameters 
determinate. The formulas for the downwash induced by a doublet distribution on super- 
sonic leading edge panels are derived in section A1 of appendix A. 

For the subsonic leading edge panels, we use a doublet distribution of the form given in 
equation (9) for the supersonic leading edge; namely, 

p = (x - e 2 y) (a + bx + cy + dxy) (10) 
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Here the origin is at the upstream corner of the panel and x - e^y — 0 defines the panel 
leading edge. The doublet distribution of equation (10) is continued along strips in the x 
direction by its value on the downstream boundary at x - and since p vanishes on the 
leading edge, p also vanishes for y > S 2 (see fig. 9). The range of integration for control 
points in the three distinct regions of influence of the panel is illustrated in figure 9. The 
four free parameters in equation (10) are determined by imposing the edge conditions and 
the downwash condition. With the choice of doublet distribution given in equation. (10), 
the condition of vanishing doublet strength on the leading edge is automatically satisfied. 

The continuity of the doublet gradient across the downstream edge is enforced when 
determining the parameters for the interior panels. The continuity of the doublet strength 
along the upstream panel edge determines two parameters (a and b) in equation (10). One 
relation is to satisfy the continuity of the normal gradient of the doublet strength on the 
downstream corner of the same edge, leaving only a single parameter to be determined by 
the downwash condition at one control point, see figure 6. In case the upstream Mach line 
edge is a special Mach line, the continuity'of the doublet normal gradient is relaxed and 
replaced by a downwash for the same reason described for the interior panels. 

The derivation of the downwash from the subsonic leading edge panel is presented in 
section A3 of appendix A and in section -B5 of appendix B and results are given by equations 
(A 13) and (B19). The calculated downwash on the subsonic leading edge has a logarithmic 
singularity as seen from the function f^ in equation (B5). Since p is approximated by a 
polynomial, the tangential velocities are finite on the subsonic leading edge whereas analytical 
solutions for subsonic leading edge wings contain an inverse square root singularity in the 
tangential velocity at the leading edge. No wakes were included in the examples computed 
by the doublet Mach line panel, since there was no subsonic trailing edge. For all trailing 
edges, continuity of p on the wing with the trailing vortex sheet must be maintained. From 
equation (7), we see that for continuity of pressure across the sheet, the doublet distribution 
on the wake must be of the form 


M = P n (x-Y) 

where P n is a polynomial. From equation (9) the expression for p on the panel adjacent 
to on either a supersonic or subsonic trailing edge of the form 


x = e 3 y 


is a polynomial of the fourth degree in x or y. Hence, P n must also be a quartic in its 
argument x-y and the coefficients must be chosen to make p continuous at the trailing 
edge. For a subsonic trailing edge panel a downwash control point must be located on the 
trailing edge to satisfy the Kutta condition of smooth flow. 

As seen from equations (8), (9), and (10) for the doublet distribution on a panel, there are 
four parameters to be determined by appropriate conditions on each panel. For supersonic 
leading edge panels we require the polynomials 


a + bx + cy + dxy 



to agree with adjacent leading edge panel .values at all comers. This ensures continuity of 
the gradient along all straight supersonic leading edges and leaves three parameters to satisfy 
the downwash boundary conditions at three control points. 

For interior panels, continuity of the normal derivative.(3/0x or 0/0y) of the doublet'strength 
is applied at the downstream ends of the two leading edges of the panel, leaving two 
parameters free to be determined by downwash boundary conditions at two control points 
on the panel. However, for interior panels in which one of the leading edges is a special 
Mach line along which discontinuities in velocity or velocity gradient may occur, the require- 
ment of the continuity of the doublet gradient is relaxed on this edge and is replaced by a 
downwash control point (see fig. 6). 

For subsonic leading edge panels, two relations are required to satisfy continuity of p- with 
the upstream Mach line edge. One relation is used to satisfy continuity of the normal gradient 
of the doublet on the downstream comer of the upstream Mach line edge, leaving only a 
single parameter to be determined by a downwash control point. 

The control point locations and continuity conditions are illustrated in figure 6 for the wing 
planform used to test the method. Design type boundary conditions also are possible with 
the xg component of the perturbation velocity prescribed at the control points instead of 
the wing camber slopes. 

4.2 DISCUSSION OF RESULTS FROM THE PLANAR 
MACH LINE DOUBLET PANEL METHOD 

To test the method described in the foregoing discussion, a cambered swept wing of zero 
thickness having a parabolic arc profile was chosen. The camber is defined by the formula 

zj = 2t (l - Cyj) (xj - x x 2 ) (1 1) 

where the y j axis is along the leading edge and the xj axis perpendicular to the y j axis. A 
closed form exact solution of the linearized differential equation was derived to check against 
the panel method for the portion of the wing planform described in figure 6. 

For the region downstream of the supersonic leading edge and upstream of the special Mach 
line, the panel method yields results in close agreement with the exact solution from 
linearized theory for C = 0.3 and a wing sweep of 30° at a Mach number -</2~as shown in 
figure 1 0. The pressure distribution normalized to maximum camber was fairly insensitive 
to location of the control points in the panels. For the supersonic leading edge triangular 
panel, best results were obtained for the interior control point located at the downstream 
comer. Increasing the panel density improved the pressure distribution. 

In the region behind the subsonic leading edge (or free edge aligned with the freestream), 
the predicted pressure distribution was poor. Increasing the panel density only served to 
increase the deviation from the closed form solution. A typical pressure distribution along 
a line close to and parallel to the special Mach line emanating downstream from the comer 
is shown in figure 1 1 . The difference between calculated downwash distribution and wing 
slope follows a similar oscillating behavior in figure 12. The zero downwash difference 
values occur at control points where the boundary conditions were applied. 
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Continuity of'the normal derivative of the doublet distribution at comers ensures the con- 
tinuity of the pressure at corners. The pressure distribution on both sides of a row of panel 
edges as shown in figure 13 demonstrates that this condition was properly applied in the 
solution. Tire greatest discontinuity in pressure appears at the middle of the panel edge. 

To test whether the difficulty stemmed from the comer where the subsonic and supersonic 
leading edges join, the flat plate conical flow solution was added to the panel method- to ' 
eliminate the singularity along the special Mach line. The pressure distributions along- 
Mach lines in the region behind the free edge are shown in figures 14 through 17. Considerable 
improvement resulted and the match with the exact solution in first column of panels down- 
stream of the special Mach line is excellent. Away from the comer, the pressure distribution 
appears to oscillate about tire exact solution. From the form of the oscillations if would 
appear that they are induced by the triangular subsonic leading edge panels, which for this 
case is a free edge. 

A few attempts were made to alleviate this difficulty. An additional control point was 
added to the triangular leading edge by adding to the doublet distribution the term 


p = cx [(y/x) - e 2 (y/x) 2 ] (12) 

which corresponds to a conical flow pressure distribution. This failed to improve the 
pressure distribution. 

Adding control points in the panels downstream of the special Mach line and applying mean 
square solution techniques to the set of simultaneous solution had a small effect in 
smoothing the result. However, as the condition of the continuity of normal derivative at 
panel corners was relaxed, the discontinuities of pressure along the panel edges become larger 
and also affect the solution in the region behind the supersonic leading edge. 


11 



5.0 PLANAR MACH LINE SOURCE PANEL METHOD WITH 
LINEARIZED BOUNDARY CONDITIONS 
5.1 THEORETICAL.DESCRIPTION OF THE METHOD 

To solve the thickness problem for a planar wing, we use a source distribution over the wing 
planform. The velocity potential in terms of Mach line coordinates of equation (3) may be 
written down by appropriately dropping the z derivative in the expression for the doublet 
potential. Thus, we have 


1 ff a(ta) dldr? 

27rM S w \/(x-U(y-T?)-M^z 2 /4 (13) 

This integralhas the property that 

h ' z = ±0 = ±P (X ’ Y)/2 = ( dz/dx o) z = ±0 

where dz/dx 0 is the wing thickness slope. Dropping the factor 2 relates a directly to the 
wing upper slope. The pressure on the wing is given by 

C p = -20 Xo = -M (^ x -0 y )/Bl z=o (14) 

Since we are interested in the flow properties on the wing, it is convenient to set z = 0 
in equation (13) before differentiating with respect to x and y. Thus on the plane z - 0, 

* = __L /' y(x) f m o&nimL 

* u / i VCx-iKy-n) 

\/y-y(x) Vx-e(y-7j 2 ) 

=- 7rE f f a ( x -$ 2 >y~n 2 )mri (15) 

o o 

where x = 2(y) or y - y(x) describes the wing (or panel) leading edge. 

The source panel method follows essentially the same procedure as the doublet panel method 
and we will use the same wing planform paneling to illustrate the method (see fig. 18). For 
supersonic leading edge panels, a quadratic source distribution in the following form is used: 

o = a + bx + cy + dxy (16) 

where the origin is at the upstream comer of the panel. One parameter is fixed by requiring 
o to be continuous along the leading edge. The other three parameters are determined from 
tangency conditions at three control points. The source strength a is continued outside 
the panel along characteristic strips by its linear values on the downstream panel edges in 
the same manner as described for the doublet strength. 
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Continuity of source strength with interior-panels is maintained by defining the source 
strength on interior panels by a relation of tKe form 

a = dxy (17) 

where the origin of coordinates is at the upstream corner of the panel. Since there is only 
one parameter d, it is determined by a single downwash control point (orpressure).at the 
panel center. The source strength of the interior panel is continued by its linear relation on 
its downstream boundaries in a manner similar to the doublet strength distribution. 

For subsonic leading edge panels, the source distribution is given by 

a = cy + dxy (18) 

where the two parameters are determined by imposing the tangency condition at two 
distinct points shown in figure 1 8. Unlike the doublet panels, the value of a on the up- 
stream characteristic strip is continued into the subsonic edge panel since a does not neces- 
sarily vanish on wing edges, either subsonic or supersonic. 

For the planform in figure 1 8, the relation for a on a given panel is derived in appendix A 
and appendix B along with formulas for the pressure coefficient. The same basic functions 
appear in the pressure coefficient for the source as in the formulas for the doublet downwash. 

Either analysis (i.e., upper surface slope) boundary conditions or design (i.e., pressure 
coefficient) boundary conditions may be applied. The location of control points is shown 
in figure 1 8. Since Cp has a logarithmic singularity on the subsonic leading edge, the control 
point must be moved off the edge of the subsonic leading edge panel for design boundary 
conditions. 

5.2 DISCUSSION OF RESULTS FROM THE PLANAR MACH LINE 
SOURCE PANEL METHOD 

To test the source panel method, a wing with a symmetric parabolic arc profile was chosen 
which has the same spanwise and chordwise variation in slope for the upper surface as the 
cambered zero thickness wing used to test the doublet method (eq. (11)). The angle of 
sweep is 30°, and the freestream Mach number is-y/2^ The subsonic leading edge angle 
was selected to make an isosceles triangle with the downstream Mach line boundary as 
shown in figure 19. The exact linearized solution is easily obtained in closed form and 
was used as a comparison with the panel method results. 

Figure 19 shows a plot of the pressure distribution from the panel method using six panels. 
The agreement with the exact solution is very close. Although the results of the panel 
method is indicated by circled points, the distribution is continuous in the whole region. 
Figure 20 shows a comparison of source distribution from the panel method with the wing 
slopes along the lines shown in the planform in the figure. This indicates that the boundary 
conditions are appropriately satisfied. 


13 - 



The design mode was also tested using the same wing. In this case the applied boundary 
condition was the theoretical pressure distribution from the exact linearized solution and the 
coefficients of the source distribution were solved for. The resulting source distribution is 
seen. to .be-in close-agreement- with- the-wing'slopes as indicated in figure 21. 
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6.0 CONCLUSION 


The source panel method yielded very accurate results when either analysis (downwash) or 
design (pressure) boundary conditions were applied. The doublet method gave accurate 
results for analysis boundary conditions on regions behind a supersonic leading edge. Both 
methods were very stable and showed little sensitivity with variation in the location, of 
control points on the panel. Refining-the paneling improved the accuracy. 

However, when the-portion of the wing behind the subsonic leading edge was included in 
the solution, the doublet method became unstable and sensitive to control point location. 
Refinement of the paneling did not improve the accuracy but increased the oscillations in 
the panel solution and its departure from the exact linearized solution. The difficulty 
appears to lie in an unsuitable treatment of triangular panels occurring along the subsonic 
leading edges. - 

Research was discontinued on this method before the difficulty on the subsonic leading 
edge was resolved. The success and stability of the subsonic method for nonplanar surfaces 
developed by Johnson and Rubbert (ref. 1) suggested that it would also be suitable for 
supersonic flow. Such an approach also has the advantage of producing methods for both 
subsonic and supersonic flow which are compatible with respect to paneling and control 
point location. 
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APPENDIX A 

DERIVATION OF BASIC FORMULAS 
FOR THE PLANAR MACH LINE PANEL METHODS 


A1 FORMULAS FOR THE DOWNWASH FROM SUPERSONIC 
LEADING EDGE DOUBLET PANELS 

We shall now develop the basic formulas for the downwash from the panels located at three 
different regions on the planform using the doublet distributions presented in section 4. 1 . 
Consider first the downwash induced by a supersonic leading edge panel. Since the doublet 
strength is continued on Mach line strips by its value on panel edges, the range of integration 
for points on and off the panel is illustrated in-figure-8. 

For simplicity we shall consider cases b and d only. The other cases can be easily obtained 
by dropping some of the terms. Thus, we have from equation (4) for case b when x < 0, 
y>sj 


w = 


m a 2 


2 ir 3xdy 


-x/ei -ejr? sj " -eisi 

f f M(£,n)dfdT?/H 0 + f f m(€,si) d?dij/H 0 

S| x y x 


(Al) 


and for case d when x > 0, y > s^ 


w =• 


M d 2 


2?r 9x3y 


f 


l s l 


r / .- € i 7 ? o 1 

[ o / ! f ) «J 


dr? 


♦ /' [ / £ ‘ S1 | iS£i ) dt + / eiSa’dfl 
0 


1 S 1 


H. 


H r 


J 


dr? | 


(A2) 


where 


h q = 


With the integrals in their present form, the integration must be performed before the 
differentiation. When the variables of integration are changed to 


, J7'=vy -t? . ( A3 ) 

the differentiation may be performed first. The integral equation (Al) then takes the 
form for case b when x < 0, y > Sj 
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a 2 

9x3y 


v/y + x/e 1 

< / 

yy- s i 


■ v/x + 6! (.y-T? 2 ) 

/ M (x-£ 2 ,y-77 2 )d$dr? 


+ s 

o 


yy - Si ^/x + ejSi 


f 


m(x~£ 2 , X! ^d£drj 


Performing the differentiation and noting that m vanishes on the leading edge; i.e., 
m(- eiy,y) - 0 leads to 


W = 2M 

it 


f 1 My (-ei (y-t? 2 ),y-t? 2 )di?/2i/ x + e 1 ( y_7j2 ) 

Vy-sj 


Vy + x/ej 

+ f 

\/y-sj 



for case b in which x < 0 and y > Sj. Similarly, equation (A 2) becomes 


w = — ^ 

7T 


Vy ' 

/ My ^y-i? 2 ^y~77 2 )d77/2Vx + e 1 ^y--rj 2 j 

vy- s i 


(A4) 


Vy" v/x + e! Jy-rrj 

+ f /M X y (x-| 2 ,y-i? 2 ) d|di7 ► 

Vy-si Vx- 
for case d in which x > 0, y > si 

The cases a and c in figure 8 easily result from the preceding equations by setting 
Vy - sj “ 0 whenever y < $i . 

Evaluating the integrals using equation (9) for ]i leads to the downwash in the form 

W = (MVej/^ jji gj(x,y) + b g 2 (x,y) + cg 3 (x,y) + d g 4 (x,y)j 
The functions gj are presented in section B3 of appendix B. 


(A5) 


(A6) 
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A2 FORMULAS FOR THE DOWNWASH FROM INTERIOR 
DOUBLET PANELS 


For the interior panel, the regions of integration are depicted in figure-4. We-shall- consider 
the-case y > sj and x > ejSj. The downwash integral then becomes 


W = 


M 9 2 
27r 9x3y 


/ 


S 1 


r 0 r 0 

J j y. (tri . 

s i e i s i 

pH n («i»i^i) ii+ 


£ 1 S I ”° 


I dr? 


+ / I 11 ** ( e l s l j7? ) d £ d7 l/H 0 

S j X 


(A7) 


In terms of the variables in equation (A3), we obtain, after performing the differentiation, 


w = ^jr J J ^x v ( x '£ 2 ,y-i7 2 ) d £dT7 (A8) 

y/y- s i V^-ejS] 
since y(o,sij = 0^ =p y ^0, sj j = 0. 

The remaining cases a, b, and c in figure 4 follow by setting Vy-sj equal to zero for y < s j 
andy^TTJs] - 0 for x < e jSj, Substituting equation (8) for ju yields the following relation 
for the downwash 


w = (m V ej/^) ja hj (x, y) + b I12 (x, y) + c h 3 (x, y) + d h 4 (x, y)J 


CA9) 


where the functions hj are derived in section B4 of appendix B. 
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A3 FORMULAS FOR THE DOWNWASH FROM SUBSONIC 
LEADING EDGE DOUBLET PANELS 


For the subsonic leading edge panels, we use a doublet distribution of the form of equation 
( 10 ); i.e.. 


ju = - C 2 y j (a + bx + cy + dxy) 


(A10) 


where the origin is at the upstream comer of the panel. The doublet distribution of equation 
(A1 0) is continued along strips in the x direction by its value on' the boundary x = € 2*2 anc * 
vanishes for y > S 2 - The special cases for the regions of integration are illustrated in 
figure 9. The most general case c yields the down wash formula 




9x3y 




(All) 


In terms of the variables in equation (A3), these integrals become, after differentiation, 

yy 


W = 2M 
7 r 


/ My (e 2 (y-v 2 ),y~ f? 2 ) <fc?/2 y/x -e 2 (y- 1? 2 ) 

[ Vy-S2 


+ 


Vy \4-e 2 (y- 1 ? 2 ) 
f f M Xy (x - i 2 > y - V 2 ) d£drj 

Vy~ s 2 V* - e 2 s 2 ■ 


' (A 12) 


The cases a and b are easily derived from equation (A12) by setting yy - S 9 equal to zero 
whenever y < S 2 and similarly y'x - ^2 e< 3 ua l to zero for x < 62 s 2 - Substituting equation 
(A 10) for m in equation (A 12) and integrating yields the following form for the downwash 


W 


= ( M j^aSj (x, y) + bS 2 (x, y) + cS 3 (x, y) +-dS 4 (x, y)J 


(A13) 


where the Sj <J = 1 to 4) functions are derived in section B5 of appendix B. 
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A4 DESCRIPTION OF PANEL AND PARAMETER NUMBERING SYSTEM 


To test the method, we choose a wing planform consisting of a straight supersonic leading 
edge followed by a straightsubsonic leading edge and consider the-region-bounded-by the 
two leading edges and a Mach line which intersects the two edges as shown in figure 22. 

The planform in figure 22 is divided into two regions by the special Mach line from the 
corner which in turn are divided'by the other family of Mach lines into rectangular and 
triangular panels. The panel numbering follows the sequence as shown in figure 22. The 
ith column has i panels for i < i c . Summing theprevious panels, we see that the panel 
number N p for the kth panel in the ith column is 

N p (i,k) = i(j - l)/2 + k,.k = 1 to i. (A14) 

Since there are four parameters (a,b,c,d) associated with each panel, the parameter 
numbers associated with the kth panel of the ith column is given by 

N(i,kj) = 2i(i - 1) + 4(k - 1) + j, j = 1 to 4. (A1 5) 

Since the doublet distribution in each panel is referred to the origin of coordinates at the 
upstream comers of each panel, we define panel variables Xj, yj by 


X{ = x + i ejsj, i= 1, 2, . . . 

(A 16) 

yj = y-jsi,j<i c ,j= 1, 2, . . . 

(A17) 


where sj is the column width behind the supersonic leading edge. 
For j > i c , we have 


yj=y-i c s 1 -(j-i c +l)s 2 (A18) 

where S 2 is the width of the subsonic Mach line columns. Since e^s^ = £ 2 S 2> ^ ie definition 
of Xj for i > i c remains the same. 

A5 THE DOUBLET DISTRIBUTION FOR PANELS IN COLUMNS 
BEHIND THE SUPERSONIC LEADING EDGE (i < i c ) 

It is convenient to define the function 

G(i, k, x, y) = X N k , i ) + x % (i, k, 2) + V X N(i, k, 3) + x ^ X N(i, k, 4) . ( A19 > 

where the X N correspond to the parameters (a,b,c,d) in equations (8) through (10) and 
N(i,kj) is defined in equation (A15). 
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Then the contribution to the doublet distribution from variables related to the kth panel in 
the ith column is, for i < i c , 

** = x i-k+l yi-1 G (i> k > yw) (A20) 

for i > k > 2 for interior panels and 

M = ( x M + e x y M ) G (i, 1, x^, y^j j (A21) 

for supersonic leading edge panels. 

We shall now construct p in the kth panel of the ith column for columns behind the 
supersonic leading edge. The upstream panels in the same row as the kth panel of the ith 
column are 


N p (i-j, k-j) for j = 1,2,3, . . ., k-1 

These panels contribute to the x variation in the form 

/ \ k-2 • 

x i-k+l G(i-k+l, 1, Xj_ k , sj) + 2 Xi _ k+1 Sl G(i-j, k-j, sj 

' 7 j=l ' 

Similarly, the upstream panel numbers in the same column are 

Np(i,k-j),j = 1,2,3, . . ., k-1 
and they contribute to the variation with y in the form 


k-1 


= 1 yj_i G(i, 1, 0, y M ) + ^ ej Sj y M c(i, j, e^, y^i) 


(A22) 


(A23) 


(A24) 


The corner values from the remaining upstream supersonic leading edge panels contribute 
a constant term of the form 


k-1 

Z 

j=2 


ejSj G(i-j + 1, 1, 0, 


(A25) 


Contribution to the constant from the interior panels upstream of the k - 1 row from 
panels numbered 


Np(i-j,k-j-l), j = 1,2,3, . . ., k-3 
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IS 


k " 3 9 / \ 

z: ejSj Gli-j, k-l-j, ejSj, Sjl 

j=! ' . 1 


Contribution to the constant from upstream k - 2 row is 

k-4 


2 e l s 1 ^ G^i-j, k-2-j, ej, sj, sjj 
j=l 


and the total contributions to the constant from all upstream interior panels are 

k-3 k-2-n 


Z Z e l s i G(i-j, k-n-j, eis^, sj ) 
n=l j=l v ' 


Combining.equations-(A23), (A24), (A25), and (A26) with the panel function itself 
yields the following value of for the kth panel-in the ith column 


Mik = x i-k+l ^i-l G (h k » x i-k+l> y i-l) + x i-k+l G ( i_k+1 > l > x i-k» s l) 


k— 2 

S 1 x i-k+l g (h, k-j, Xj_ k+1 , sj) + e iyi _f c(i, 1, 0, y M ) 


+ 2 e l s l Yi-| G (i* j> e l s i > Yi-l) +2 e l s l G (i-j+l. 1,0, sj) 
j=2 v ' j*2 


k-l 


k-3 k-2~p 

+ 2 2 • € i s i Z Gli-j, k-n-j, ejs^ sj 

n=l j=l v 


(A26) 


(A27) 
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Equation (A27) holds for k > 2, if those summations are dropped for the particular values 
of k in which the upper limits are less than the lower limits. 

A6 DOUBLET DISTRIBUTION FOR PANELS IN COLUMNS 
BEHIND THE SUBSONIC LEADING EDGE (i > i c ) 

Let i c be the last column for which the leading edge is supersonic. The total number of 
panels in the region to the left of the special Mach line is 


icOc" 1 - 1 )/ 2 

The number of panels for i > i c in the ith column is 

i c-( i - i c"0 = 2i c- i+1 

By summing the panels in the columns up to the ith column, we find that the panel number 
for the kth panel in the ith column is 

N p (i,k) =|^i c (i c + l) + (3i c -i + 2)(i-i c -l)j /2 + k (A28) 

The numbers corresponding to the parameters in the panel are 

N(i,kj) = 4N p (i,k) - 4 + j, j = 1 ,2,3,4. (A29) 

For the subsonic leading edge panel we have 

K il = ( x 2i c -i+l - e 2 yi-l) G (u, x 2 i c -i+l> yi-l) (A30) 

The doublet distribution cannot be made continuous by analytically continuing the 
distribution from the upstream panel downstream into the subsonic leading edge panel, 
since ju must vanish on the leading edge. Thus ju must be matched in the subsonic leading 
edge panel on the upstream edge by two relations'obtained by equating the -two coefficients 
of the x variable. 

Consider the kth panel in the ith column for i > i c . Then the panel numbers for the 
panels in the same row for i > i c are 

Np(i - j, k + j), j ~ 1 ,2,3 , . . i - i c - 1 

The remaining interior panels upstream of the special Mach line in the same row are 
N p (i c + 1 - j, k + i - i c - j), for j = 1, 2, 3, 4, . . k + i - i c - 2 

where Np(i,k) = i(i-l)/2 + k for i < i c . 
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The supersonic leading edge panel in the same row as the kth panel in the ith row is 

N p (2i c -i~k+2,l) 

Contributions to the x variable alone in fi for the kth panel of-'the ith column are given by 


i-i c -l 

X 2i c -i+2-k s 2 G ( i_ h k+j, x 2i c -i+2-k’ s 2) 
k+l-i c -2 

+ 2 x 2L-i+2-k S 1 G ( i c + H> k+i-i c -j, x 2i J+2-k, S 1 

j=l c X C 

+ x 2i c -i+2-k G ( 2i c -i_ k +2 > 1> x 2i c -i-k+l> s l) 


(A31) 


Contributions to-the y variable in p for the kth panel of the ith column are 


k-2 

Z 

j=l 


e 2 s 2 ^i-l G (h k-j, e 2 s 2 , Yi_i) 
" e 2 y i G (i’ 1> e 2 s 2> yi-l) ’ 


+ x 2i c -i+2-k y i-l G (h k, x 2ic _i +2 _ kj ‘ (A32) 

Except for a constant, is given by the sum of equations (A3 1) and (A32). For k = 2, 
this constant will be zero since vanishes for y:_j “ s 2 and x 2 j = e 2 s 2 , a point on the 
leading edge. > ( ' 

To make n continuous on the upstream edge of the first panel in the ith column, we set 
k = 2 and replace i by i-1 in equations (A3 1) and (A3 2). 


Combining the results for yj_j = 0, and equating to p defined by equation (A30) for yj_j = 0, 
we obtain 
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x 2i c -i+l G ( i > 1. x 2i c -i+l’°) 
i— i c — 2 

” 2 x 2i„-i+l s 2 G (H- 1 » 2 ^‘> x 2L-i+l » *2) 

j=0 c ' c . / 

i-i c -1 . . 

+ * x 2i c -i+l S 1 G ( i c +1 -j> M c + H, x 2i c -i+l’ s l) 

+ x 2i c -i+l G ( 2i c“ i+1 * x 2i c -i, Sl ) (A33 > 

For the corner panel i = i c + 1, equation (A33), takes the form 



Using the relation x-i; _j_ 1 = Xo.- -ei si , combining coefficients of Xy. and x 2 2ir-i> and 
setting the coefficients in equation (A33) equal to zero yields two relations which reduce 
the free parameters of the subsonic leading edge panel to two. 

To determine the constants resulting from the upstream panel comers for the columns 
i > i c , we note that the constants from the panel comers are incorporated into the subsonic 
leading edge by means of the continuity conditions in equation (A3 3). Hence, for the kth 
panel in the ith column, the contributions to the constant term in occur only for the 
rows k - 1 to 2. Since the constant term itself does not affect the pressure distribution it 
will not be written down. Accordingly, the doublet distribution for i > i c takes the form 

**ik = x 2i c -i+2-k y i-l G ( i > k » x 2i c -i+2-k, y i-l) 

^c- 1 . 

+ 2 j x 2i c -i+2-k s 2 G (i-i> k ^> x 2i c -i+2-k, s 2j 


k+i-i c -2 

+ 2 1 x 2 i c -i+ 2 -k S 1 G (i C +1 -j’ k+i-i c -j, x 2 i c -i+2-k’ s l) 

+ x 2i c -i+2-k G (2i c -i-k+2, 1, ^2i c -i-k+l’ s l) 

k ~ 2 / \ 

+ 2 e 2 s 2 y i-l G i > e 2 s 2 , Yi-i) 
j=l ' ' 

- €2 y j-2 G (i. 1 , e2 s 2» Yi-l) + constant 


(A34) 
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A7 FORMULAS FOR THE DOWNWASH IN COLUMNS 
BEHIND SUPERSONIC LEADING EDGE (i < i c ) 

Let be the contribution to the downwash b<j>/bz from .the.panel variables-in-the-kth 
panelofthe ith column. For i < i c , is defined in equation (A6) for a supersonic leading 
edge panel (k = 1), in equation-(A9) for interior panels and by equation (A13) for subsonic 
leading edge panels. For the total downwash at the point on the kth panel of the ith 
column we add the contributions from all the panels upstream of the kth panel in' the ith 
column. The contributions from the same row are 

k-1 

2 W i-m,k-m 
m=0 


The contributions from the row above the kth panel are 


k-2 


£•' '^i-m.k-m-l • 

m-0 


and finally from all panels upstream .of the kth panel of the ith column we have 

-k-1 k-n-1 

W = 2 2 - ^i-m, k-m-n 

n=0 m=0 

Since the downwash W is computed from a different formula for the supersonic leading 
edge panels than for the interior panels it is convenient to separate them, thus we write 


k-1 
W= 2 W; 


k-2 k-n-2 
+ 2 


i-m,l" rz ' ^ ^i-m, k-m-n 

m=0 . - n=0- m=0 

Substituting equations (A6) and (A9) into equation (A3 5) yields 

■4 k-1 ' 

5rW/M N /?= 2 2 X N(i-m,l j) g j ( x i-m-l 5 y i-m~l) 

j=l m=0 v > 

X N(i-m, k-m-n, j) h j ( x i-k+n+l> y i-m-l) 


4 k-2 k-n-2 

+ 22 2 
j=l n=0 m=0 


(A3 5) 


(A36) 


where N(i,kj) = 2i(i-l) + 4k + j - 4. 


26 



A8 FORMULAS FOR THE DOWNWASH IN COLUMNS 
BEHIND SUBSONIC LEADING EDGE (i > i c ) 


We consider the downwash in the kth panel of the ith column where now i > i c . The panel 
number in the i c column in the same row as the Np(i,k) panel is 

N p (i c ,H c +k-l) 

The contribution to W for all the panels in columns behind the supersonic leading edge 
is found from equation (A3 5) by setting i - i c and replacing k by i-i c + k-1 or . 


i-i c +k-2 


r 

m=0 



i-i c +k~3 
+ Z 
n=0 


i-i c +k-3-n 

^ . ^ L-m ,i-i r +k- 1 -m-n 

m=0 c . c 


The contributions of the downwash to the kth panel of the ith column from panels in the 
columns whose numbers are greater than i c and in the same or k row are 

i- i,, - 1 

2 ^i-m.k+m 
m=0 


Similarly for the k-1 row 

i-ic-i 

2 ^i-m.k+m-l 
ra =0 

and for all k-2 rows 


k-2 i— i c — 1 

2 2 \y. , 

n=0 m=0 w i -m > k+m-n 


The contributions from the subsonic leading edge panels are 


i-i c -! 


2 

m=0 


Wi- m ,l 


■and from the remaining interior panels 


i— i c — 2 i-i c -2-n 

2 2 ‘ W i-m-n-l ,2+m 

n=0 m=0 
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The total downwash fora point on the kth panel of the ith column for-i > i c is, finally, 


i-i c +k-2 M c +k-3 

W = 2 + * 

m=0 n-0 


i-i c +k-3-n 
2 

m=0 


■^i c -m,i-i c +k-l -m-n 


i— i_— 1 
+ 2 W 
m=0 


i-i c -2 
+ 2 
n=0 


k-2 i— i c — 1 

i-m,l + ^ ^ ' ^i-m,k+m-n 

n=0 m=0 

•i-ic-2-n . 

^ ^i-m-n-l ,2+m 


For the supersonic leading edge we have 

* W i c -m,l/ M V^T 

4 

j?j -^-N(i c -m,l j)%‘ ( x i c -m-l 5 y i c -m-l) 

for the subsonic -leading edge 


" W i- 

4 

= j -1 X N(i-m,lj) S j ( x 2i c -i+m-l, y i-pi-l| 
and for the interior panels in columns i > i c 

4 

*W ijk /M/i^ = 2^ x N(i,kj) bj ( x 2i c -i-k+2, ^i-l) 


(A37) 


(A38) 


(A39) 


(A40) 
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and for interior panels in columns i < i c 

4 , 

irW ik /MVq = _^ XNakj) 1 ^ ( x i-k+l»yi-l) < A41 > 

Note that for i < i c , 

yj_l =y-(i-l)s 1 

N(i,kj) = 2i(i - 1) + 4(k - 1) + j (A42) 


while for i > i c 


Yi -1 = y- i c s l"( i c -i+2 ) s 2 

N(i,kj) = 2i c (i c + l) + 2(3i c - i + 2) (i - i c - l)+ 4(k-l) + j (A43) 

Continuity of downwash along the supersonic leading edge at panel comers requires 
continuity of the derivatives of p at panel comers. For the point between the i and i+1 
columns, this becomes 

X N(i+l,l,l) = X N(i,l,l) ~ e l s l X N(i,l,2) + S 1 X N(i,l,3) * e l s l 2 X N(i,l,4) ( A44 > 


A9 MATRIX EQUATIONS FOR THE DOUBLET PANEL METHOD 

At each panel both equations (A36) and (A37) for the downwash reduce to an equation of 
the form 


N 

2 AnXn-W 

n=l 


where W is the local camber slope dz/dxQ for the wing at a selected control point of the 
panel. Four such equations are needed for each panel since there are four variables to be 
defined at each panel. For subsonic leading panels two of these equations are.provided by 
the condition that p be continuous with the upstream panels by matching the coefficients 
of x and x^ in equation (A33). Another condition is the continuity of the normal gradient 
of p at a point on the upstream boundary of the panel. This condition is relaxed for panels 
bordering a special Mach line where pressure or pressure gradients may be discontinuous. 

The designation of downwash control points is indicated by dots and continuity of normal 
derivatives of p at panel boundaries is indicated by arrows in figure.6 with the arrow heads 
lying.in the panels in which the conditions apply. These continuity conditions are 
prescribed at the supersonic leading edge panel corners at 

x = 0, y = si 


29 



for interior panels at 


x = 0, y = s l 
x = e 1 s 1 ,y = 0 

and for subsonic leading edge panels at 

x = e l s l =e 2 S2,y = Q, 

‘when expressed in local panel'co ordinates. Noting that 


and 


[ x 2y m G 

Xg,y m j 

& 2x £> y m ) 


[ x ey m G l 

[*> x £>y m )] 

= x £ G (oc, p, Xg,2y m ) 



(A45) 


the conditions that 3p/3x is continuous at x^^j = ejSj and y^_j = 0 for the kth panel in 
the ith column is ' • ‘ • 


_ 3 _ 

3y M i,k 


3^ ^i-l,k-l 


ejSj, sj 


e l s l’0 

4 - 

From equation (A27) with the aid of equation (A45), equation (A46) becomes 
sj G^i,k, e2S|,0^+G(i,'l',0,0) + Sj 2 G^i, j, e.js^oj 

= s i G^-l,k-l,e 1 s i ,2s 1 ) + g({- 1, 1,0; 2 sj) ’ 

k " 2 /. .. \ 

+ $j2 Gl i-1, j, e2 s l> ^ sjl 

j~2 v ' 

This leads to homogeneous equations of the form 

‘ ' / * 

N 


(A46) 


(A47) 


2 A n X n = 0 
n=l 


The other' continuity condition 


3x ' M i»k I • 3x M i,k-1 

0, S! 


(A48) 


e l s l’ s l 
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with the aid of equations (A27) and (A45) leads to 

G (i, k, 0, sj) + X N(i _ k+1> l} d - ejsj X N(i _ k+lj 1} 2 ) 

+ s 1 [^ x N(i-k+l J 1, 3)“ e l s l x N(i-k+l, 1,4)] +S! ^ G (i-j, k-j, 0, 

= Sl G(i, k-1, 2 ei s h + X N(i _ k+2} l, 1) + *l s l x N(i-k+2, 1, 2) 

+ S 1 [ X N(i -k+2, 1, 3) + e l s l X N(i-k+2, 1, 4)] 
k - 3 /. . 

•+si I G (i-j, k-l-j, 2 e^sj, sj 
j=l V 


The case for which k = 2 , must be treated separately. Equation (A 46 ) yields 

^ [ x m yj-i G (h 2 > x i-i> y-i-i) + e i Vi-i G (‘> l > °> yi-i)] e i s i > 0 

= 1 1 [ ( X i -2 + e l y l- 2 ) G ( M ’ l > x i- 2 > y i— 2) ] o, Sj 
or 

Sj G(i, 2 ,e lSl ,o) +G(i, 1 , 0 , o) = G(i-l, 1 , 0,2 8 ^ 

Similarly, equation (A48) yields 

ST [ X M y M G ( i5 2 > x i-l’ y i-l) + x i-l G (i- l > l > X i- 2 > S l)] Q g 

= [ ( x i-l + e l yi-l) G (** Uxj_ x ,yi_i)] ^ 

or 

sj G(i, 2 , 0 , sj^ + G^i- 1 , 1 -ejSj, sj^ 

= j^X N(i> 1} 2 ) + S 1 x N(i,.l, 4 )] + G (h °> s l) 


(A49) 

(A50) 

(A51) 

(A52) 

(A53) 
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Continuity of 3/a/3y on the upstream boundary of the subsonic leading edge panel of the ith 
column is found by setting k = 2, replacing i.by i - 1 in equation (A34), and differentiating. 
Thus with equation (A30), we have 

4r[ x 2i c -i+l y i-2 G ( i_1 > 2 » x 2i c -i-l’ y i-2) 

- e 2 yj_3 1. e 2 s 2> y i-2) 1 

' Je 2 s 2’ s 2 

= 3 7[( x 2i c -i+l ' e 2 y i-l) G 0> x 2i c -i+l’ y i-l) ] ^ Q 

or ' ' 

s 2 G (i - 1> e 2 s 2 , 2 s 2 1, e 2 s 2 , 2 s 2 ^ 

= s 2 [ X N(i,l,3) + e 2 s 2 X N(i, 1,4)]" G ( i ’ l ’ e 2 s 2> °) 

Equations corresponding to equations (A47) and (A49)-for interior panels for i > i c can be 
derived in a similar manner using equation (A34) for 

A10 SOURCE DISTRIBUTION FOR COLUMNS BEHIND 
• ' ; SUPERSONIC LEADING EDGE 

The source panel method follows essentially the same procedure as the doublet method, and 
we use the same wing planform to test the method. For supersonic leading edge panels, a 
quadratic source distribution in the following form is used: 

1 a = a + bx + cy + dxy '• (A54) 

where the origin is at the upstream comer of the panel. One parameter is fixed by requiring 
a to be continuous along the leading edge. For the supersonic leading edge panels denoted 
by subscripts 1 and 2 and numbered 1 and 2 in figure 22, the continuity- conditions take 
the form 


a2 - a i -bei s i + c$i - deiSi 2 (A55) 

for a leading edge of the form x + e jy= 0. Here sj is the panel width (y direction). 

For interior panels, we associate a single free parameter and modify the source distribution 
by 


a = d xy (A56) 

where the origin of x and y is at the upstream leading edge comer of the panel. Since it 
vanishes on the leading edges, it preserves continuity when a is continued downstream along 
Mach line strips by the value on the panel trailing edge in the same manner as the doublet 
panel method. 
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To construct the value of the source at some kth interior panel in column i we consider 
the same wing as in the doublet method. From figure 22 we define a in the leading edge 
panels as 


a l ~ a l + b l + c l y + dj xy ’ 
a 2 ~ a 2 + b 2 X 1 + c 2 y ! + d 2 x iyi 

04 = a 4 + b 4 x 2 + c 4 y 2 + d 4 x 2 y 2 (A57) 

where Xj,yj has the same meaning as in the doublet paneling. To preserve continuity the 
• source distribution in panel 3 takes the form 

°3 = a i( x ’ s l) +a 2( 0 ’ y l)“ ff 2( 0 ’°) + d 3 x l y l (A58) 

Since a 2 = a$ - b ejSj + csj - d ejSj^ from equation (A55) we have 

°T( x ’ s l) -a 2^ 0>0 ^ = ( b l +d l s l) ( x + e l s l) 

= (bi+di$i) X! 

from which 

a 3 ~ ( b l +d l s l) X 1 +a 2 + c 2 y l +d 3 X 1 y l (A59) 

and similarly 

a 5 = (b 2 + d 2 si) x 2 + a 4 + c 4 y 2 + d5 x 2y 2 (A60> 

To find ug we write 

°6 = °3 ( x l 5 s l) + CT 5 ( e l 8 l» y 2) “ a 5 («1 8 1»0) + d 6 xy 2 (A61) 

or 

°6 = ( b l + d l s l) X 1 + d 3 S 1 X 1 + d 5 e l s l. y 2 + d 6 + y 2 

+ a 4 + c 4 y 2 + (b 2 + d 2 sj) ejsj (A62) 

We note that the supersonic leading edge panels are continued in the y characteristic 
strip in the form 

(bj+djs^x! 

and the constant term, a, only enters for the leading edge panel of the same column as the 
panel under consideration. 
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The panel numbering system for the source panel method is identical to that used in the 
doublet panel method. The panel number for i < i c for the kth panel in the ith column is 
then given by 


Np = i(i— 1 )/2 + k , k = 1, 2, . . ,,i 

For k=l, the leading edge panel has four associated parameters while for all other k there 
is only one associated parameter. The variable (or parameter) number in the kth panel 
(k > 1) is 


N(i,k) = i(i-l)/2 + k + 3i = i(i + 5)/2 + k, k > 1 (A63) 

For k = 1, the variable numbers are. 

Nj (i,j)“i(i + 5)/2+j-3 , j=l, 2,3,4 

The pattern in equation (A62) can be generalized and for the kth panel in the ith column we 
have 


ff ik".( x N 1 (i-k+1,2) + S 1 X N 1 (i-k+l,4)) x i-k+l 


k-2 

+ si 2 x N(i-j,k-j) x i-k+l 
J=1 


k-2 

+ ejsj I X N(i,k-j) y i-l 
J=1 

v 

+ x N 1 (i,l) + X N 1 (J,3) y i-l + x N(i,k) x i-k+l y i-l 


k-2 
+ 1 
j=l 




e l s l 


+ 


e l s l 2 


k-3 

I 

n=l 


k-2-n 

X X N(i-n,k-m-n) 
m=l 


(A64) 
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where N(i,k) denotes tire variable number for the kth panel defined for k > 1 in equation 
(A63), and Nj(ij) denotes the variable numbers for the first panel of the ith column given 
by the equation following equation (A63). 

A1 1 PRESSURE COEFFICIENT FROM THE SUPERSONIC LEADING EDGE 
PANELS AND FROM THE INTERIOR PANELS 

The velocity potential resulting from a leading edge panel for points in the panel is seen from 
equation (15) to be 


Vy +x / e i \/ x+e i (y - v 2 ) 

= I f a ( x “f 2 .y“»? 2 )d|di? (A65) 

0 0 ' 

where x + e^y = 0 defines the leading edge. 

For points outside the panel, x < 0, y > sj, we have 


Vy+x/ej v&?) 

^ = “irM j / f oU~$ 2 ,y^l 2 \d$<in 

V y-sj 0 


+ I 

0 


V x+e jSj 

f o (x-£ 2 sj j did n 


(A 66) 


For x > 0, y < s j, we obtain 



I 




x + ei 


0 


V3T' 


7 


(y - £ 2 ) 

x-£ 2 , y-i? 2 j d|d7} 


VT 

VF 

+ / 

I o/0,y-J? 2 ) d|dr} ► 

0 

0 ' 


(A67) 
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Similarly;.fpr, x.>' 0, y‘> S] , 


'/iiJLj • y x+? i( y "f) 2 •• : x - 

' ttM: j - J. '' * J ’ y-f? 2 jd$dr? 

Vy~ s i V* ' 7 


VT r'f*—, 

+ J J o (o,y-r? 2 jd£di7 


' y ~ s i 


! i 


' »« J 


V^y^s] Vx + ejsj . • v 

+ J J a(x-£, sjJd^dT? 

0 Vx~ 


vr-si y5T 7 

+ y J g (0-,s 1 )dSd7, 
o o ' ; ■ 




(A68) 


Only equations (A66) and (A68) need be considered,, since the other two cases result from 
setting >/y - S j - 0 when y < s^ andy^ = 0 when x < 0. From equation (A66) the pressure 
coefficient becomes r 


C P “ B (*x + ^y) - " 

-V J 

(«♦«.) / 

a (r e i (y-’j 2 )>y-’? 2 )d77 


Vy"- s i 

2Vx 2 + ei ^y ‘ • 


yjy + x/ei Vx + ej (y -ij2) 

+ /' . / .. ’ px (x-$ 2 ,y-?? 2 ) + a y (x-£ 2 , y-t? 2 y)l d$d7? 

Vy-si o ■ ‘ ’ L ' 

Vy -si V x + e i s i . - 

+ / / a x (x-i 2 ;sfV.d^ > 

0 • 0 

yy- s i 

where the term J n (~ e l s l> s l)dn- * s dropped because of continuity of a at panel corners. 
® 2-v/x + e^i 



Equation (A68) yields 


c =-4- 

C P. irB 


i 


/ 1+c \ y^ ~ q (- e i (y-T? 2 )>y- 
1 0 v^T ' 


%/y - Vx + ej (y-7? 2 ) n 

+ / / [a x (x-^ 2 ,y-i? 2 ) + a y (x-? 2 > y-T ? 2 )] 


y-si >/F 


d£dr? 


Vy" Vx" ' 

+ . f f a y (0,y-Tj 2 Jd$dr? 

vy- s i 0 


Vy-si VxTqsf 
+ f f ff x (x - £ 2 , Sj ^ d£dr? 

0 \/x~ 

where we have neglected the terms 


(A69) 


1 

2V7 


VF 

f a (0,0) d£ and 
0 


1 


2-v/x + ejSj 


Vy-n 

f a (” e l s l> s l) di? 
0 ■ ' 


which are cancelled because of the continuity of a at panel comers. 
Choosing a(x,y) = a + bx + cy + dxy and integrating yields Cp in the form 

c p = (Vq/irB) [agj (x,y) + b g 2 (x,y) + c g 3 (x,y) + d g 4 (x,y)] 
where the g functions are derived in section B6 of appendix B. ‘ 

For interior panels, we use a source distribution of the form 

a - d x y 


(A70) 


(A71) 


with the origin of the coordinates at the upstream comer. The velocity potential for interior 
panels is given by ’ • 

, s/y VF 

0“-^ f f ° ( x “* 2 > y-ifid&Q (A72) 

0 0 
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for x < ejSj, y <sj. For x > ejSj and y > Sj, we have 


<t> = - 


_4_ 



v/3T . Vf V x “ e } s I 

J a {x-f 2 ,y-77 2 )'dfd7? + J J a ■) d£dij 

%/x^ei s i v^s ]' 0 


y^*i 

+ / 

0 


/ 


V3T 


e l s l 
a (ejSj, 


5j) d£d7? 


yr 

+ J J a(x-t 2 ,s 1 ) d{di? 

0 v /x - e l s l 


(A73) 


The pressure coefficient from equation (A73) becomes 


C «-± 

T> 7tB 


yr y* ■. 

/ J [o y (x-£ 2 ,y-i? 2 ) + a y (x-£ 2 , y-? 7 2 )] d^fy 

VV-si V x ~ e l s l 


v7 

+ J J o y feisj^y-r?^ d|d7? 

vy-si 0 


Vy- S i >/x 

y - J ^x'( x "^ 2 ’ s l) dfdr? 

o V x_e i s i . 


where we have used the condition a(0,y) = o(x,0) = 0. Integration yields 

C p = (4/rrB) d h (x,y) 

where h(x,y) is given in section B7 of appendix B. 


(A74) 


(A75) 
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A12 PRESSURE COEFFICIENT FROM THE SUBSONIC LEADING EDGE PANELS 

The contribution to the source distribution from subsonic leading edge panels is simpler than 
that from supersonic leading edges. With the origin of the panel coordinates at the upstream 
corner, the source distribution from parameters associated with the panel is given by 

a = cy + dxy (A76) 


Unlike the doublet panels, the value of a on the upstream panel edge is continued into the 
subsonic leading edge panel since a need not vanish on the subsonic leading edge as the 
doublet must. For points in the panel, x < e2 s 2» V s 2> we have 


•s/y* \A-e2 (y-?? 2 ) 

f I °(*-?>y-n 2 )wv 
o o 


For points outside the panel in which x > e 2$2 and y < S 2 , we have 

f 

0 


* “irM f 


Vx-e 2 (y -rr) 

f a (x-£ 2 ,y-7? 2 )d£ 

V x-e 2 s 2 


Y x “ e 2 s 2 

+ J ff^e 2 S2,y-i? 2 )d^ 


Similarly, for x > e 2 S 2 and y > S 2 , we obtain 


d7j 


* = 'iiSn / 



(A77) 


(A78) 


Vx-'e 2 (y-i? 2 ) >Jx - e 2 s 2 

f ^(x-f 2 ,y-i 7 2 )d$ + y^ a (e 2 s 2 , y - n 2 .) d{ 
V x - e 2 s l 0 


d?? 


• x/y-s 2 ^x-e 2 (y-V 2 ) 

+ J J a ( e 2 s 2 5 s 2) d $ d7 l 


0 


(A79) 
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With <p defined by equation (A79) the pressure coefficient becomes 



Substituting equation (A76) for a into equations (A80) yields a relation of the form 


c p = [cSi (x,y) + dS 2 (x,y)J (A81) 

where Sj andSj are derived in section B8 of appendix B. 

A13 CORRECTION TO THE PRESSURE COEFFICIENT FROM SUPERSONIC 
LEADING EDGE PANELS AND FROM INTERIOR PANELS WHEN THE SUBSONIC 
LEADING EDGE PANEL OF THE SAME ROW IS IN THE ZONE OF DEPENDENCE 

OF THE CONTROL POINT 

The formulas developed in sections All and A12 for supersonic leading edge and interior 
rectangular panels were derived on the assumption that the characteristic strip in which they 
lie extended all the way to the Mach lines through the control point x,y as indicated in 
figures 4 and 8. These formulas must be corrected when the Mach line cuts the subsonic 
leading edge at the end of the row as shown in figures 23 and 24. There will be an additional 
contribution from the line integral of a along the subsonic leading edge and an area integral 
subtracted from the results in sections A1 1 and A12 because of the termination of the 
characteristic strip by the subsonic leading edge panel. 

We consider a point (x,y) as shown in figure 23 and compute the pressure resulting from 
the contribution of the supersonic leading edge panel. The value of a on the downstream 
side of the triangular panel on the supersonic edge is carried on in the strip downstream of 
the panel edge by its value on the panel edge, or, as shown in section A10, 



where the origin of xj is on the upstream edge of the subsonic leading edge panel. 
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To find the correction to the pressure coefficient for the panel due to the subsonic leading 
edge in the zone of dependence of the point, we need to subtract the contribution from the 
area bounded by the lines £ -\/x7 V ~y/yl an d the subsonic edge. This integral takes the 
form 

y/y \/T 

A v'SH b+dsl ) / 

0 vx-e 2 (y-V 2 ) 


and the contribution from the integral along the segment of subsonic leading edge is 

Vy~ 

AC Pe = i! ( b + ds i) / _j2_(i-€7)(y-v 2 )dn 


- =Jl 
"'Pg irB 


x-e 2 y-n" 


(A82) 


for points y < s 2 . For y > s 2 we have, instead. 


AC P a “‘irB 


f V> 

(b + dsi) / 

_Vv~ I 


x-e 2 (y-T? 


J J dfdT ? (A83) 
0 V* - e 2 s 2 _ 


AC Pg = ?l (b + dsj) j 


e 2 ( y_1?2 ) 


• Vy- s 

■ + s 2 y . . dt ? — 
o v^x-e 2 (y-v^J 


(A84) 


The evaluation of these integrals yields a relation of the form 

AC P = jri ( b + ds l).V«2' C(x,y) (A85) 

where C(x,y) is derived in section B9 of appendix B. Note that the origin of coordinates 
for C(x,y) is that for the subsonic leading edge panel on the same row as the supersonic 
leading edge panel being considered. The correction to the pressure coefficient for rectangu- 
lar panels is found from equation (A85) with b set equal to zero. 

There is required a correction due to the constant value of a which is obtained from the 
first subsonic leading edge panel. In this panel the value of a contributed from the adjacent 
supersonic leading edge panel is 

a - a + csj + dsjx + bx 

= a-b +csj -d ejSj 2 + b+dsj x + ejSj 

This is seen to differ from that in figure 23 by the terms 

a - b ejsj + csj -d ejSj 
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The correction to Cp from these terms from the last supersonic leading edge panel for points 
downstream of the special Mach line from the corner is given by 

AC P = ir! ( a_be l s l +cs ! “ de l s l 2 ) 

“i| V^2 ^^f4( x >y) (a-bejsj+csj-dejsj^ (A86) 

for a point in the subsonic leading edge region. The function {4 (x,y) is given in equation 
(B5) of appendix B. 

We now consider the subsonic leading edge correction to the pressure from panels. In 
figure 24 we see that the contribution from the integral along the subsonic leading edge 
has the same form as o for the supersonic leading edge panels except that b = 0 and hence 
the correction to C p contains the same function C(x,y). 

A14 SOURCE DISTRIBUTION FOR PANELS IN COLUMNS BEHIND 
THE SUBSONIC LEADING EDGE 

To define the source distribution in the kth panel for the ith column where i > i c , i c being 
the last column with a supersonic leading edge panel, we consider the variable numbering for 
columns adjoining the subsonic leading edge. The subsonic leading edge panel has two • 
variables while each interior panel has only one. In numbering the variables, we start with 
the total number of variables in the columns adjoining the supersonic leading edge; i.e., 

i c( i c + 5 )/ 2 + i c 

The variable in the kth panel of the i c + 1 column is then seen to be 
N^i c +l,k^ = i c (i c + 5^/2 + i c + k + 1 for k > 1 

Each successive column (i increasing) has one less panel and hence one less variable. For 
the ith column and kth panel, we have for the variable number 

N(i,k)=? i c (i c +5^/2 + i c + 5 ; ° /i c + 2-n\ +k+ 1 

n -1 V ' 

Forming the summation and simplifying yields 

N(i,k) = i c (i c + s)/2 + (3i c ^ i + 5 ) (i- i c )/2 + k-1 

for k > 1 . 

The variable numbers for the leading edge panels are 

Nj (i j) = i c (i c + 5')/2 + (3i c - i+s) (i - i c )/2 -2+j 

j- 1,2 
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We now find by combining the contributions from all upstream columns and rows For 
the kth pa: 
column is 


the kth panel iri the ith column, the panel number in the same row and the i c or i c + 1 


i - i c + k - 1 

The contribution to the x variable from panels in the columns behind the subsonic leading 
edge is seen to be 


i-i c -! 


s 2 X N(i-j,k+j) x 2i c -i-k+2 = 2 1 x 2i c -i-k+2 

Contribution to the y variable in the kth panel of the ith column is 

' k ” 2 / \ 

• e 2 s 2 £ i x N(i,k-j) y i-l + ^ x Nj (i, 1 ) + e 2 s 2 x Nj (i,2)j ^i-l “ 2 2 Yi_1 

The constant terms resulting from the comer values of a for the rectangular panels in 
subsonic columns is 

„ i-L-2 n 
o c 

= ** V 


(A87) 


(A88) 


2 v / v 

x 3 = € 2 S 2 x x X N( i r +j,n-j+2) 

n =l j=l ' / 

k -i iV 1 2 / \ 

+ n=l f =1 ^ X N (i c +j ,i-i c +k-n-j ) 


(A89) 


The contribution to the constant from the subsonic leading edge panels is 
„ _ i-t-1 1 

S/i = 


4 52 Jl' e 2 S 2 x N 1 (i c +j,l) +x N 1 (i c U2)_ 


(A90) 


The contribution to the x variable and to the constant from the supersonic leading edge 
columns is found by setting i = i c in equation (A64) and replacing kbyi-i c + k- I and then 
substituting yj_ j. = sj . This yields 


( X Nj (2i c -i-k+2,2) + S 1 x Nj (2i c -i-k+2,4)) x 2i c -i-k+2 


i-i c +k-3 


+ si I x N(i c -j,i-i c +k-l-j) x 2i c -i-k+2 


j=l- 


i-L+k-3 


+ e lSl 2 S i x N(i c ,i-i c +k-l-j) + X Nj(i c ,l) +X 'N l (i c ,3) s l 


i-h+k-3 


+ X N(i c ,i-i c +k-l) sj x 2i c -i-k+2 + ij* ( X Nj + S 1 x Nj (i c -j,4)) e l s l 


i-i c +k-4 i-i c +k-3-^7 
+ e l s l S Z 

h=l • m=l 

s 2 5 x 2i c -i-k+2 + s 6 


Xi 


N (i c -rj,i-i c +k-l-m-n) 


(A91) 
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The source distribution oj^ for the kth panel in therith column is found by combining 
equations (A87) and (A91) or 


a ik - ( ^1 + S5 + x N(i,k) yj-1 ) x 2i c -i-k+2 
+ Z 2 yi-1 + £ 3 + z 4 + s '6 


(A92) 


Equation (A92) can be interpreted to hold for all i,k, and i > i c if those summations are 
dropped for which the upper limits are less than the lower limit or terms in which negative 
column or panel numbers occur. One exception is the first panel in the i c + 1 column 
which must be treated separately. Thus, 

%+l,l = X N 1 (i c ,l)- e l s l ?^ 1 (I C ,2) + 8 1 “ e l s l 2 x Nj(i c ,4) 


+ [ X Ni(i c » 2 ) + S l X N 1 (v4)J x i c 
+ (i c +i,i)n c + x N!(i c +i yi c 

Evaluating equation (A92) for i = i c + 1 and k = 1 .reveals that we must add the term 

. - e l*l[^ x N 1 (i c +l,2) + 8X N 1 (i xc +l,4) 

to equation (A92) for this particular panel. 

For analysis boundary conditions in which the wing shape is given and the pressure is to be 
found, either equation (A92) or (A64) is set equal'to the slope of the wing at the control 
point of the panel. These equations together with the homogeneous equations involving 
continuity of a at panel comers on the supersonic leading edge yield a set of equations to 
be solved for the x n variables. 


" 


(A93) 
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A15 COMPLETE PRESSURE COEFFICIENT FROM THE SOURCE PANELS 


The contribution to the pressure coefficient from each panel is summed in the same way as 
downwash is computed from the doublet panels. From equation (A35), we write for the 
pressure in the kth panel of the ith column for i < i c 


k-1 


P = r 


k-2 k-ij-2 


=n C i-m,l + n 5 0 2 =Q c i-m,k-m-n 


m=0 


(A94) 


where the first subscript denotes the column and the second subscript the panel number in 
the column. For the first sum, we use equation (A70) and obtain 


V^T 4 

I X 


i-m,l ttB X Nj(i-mj) g j ( x i-m-l, ^i-m-l) 


(A95) 


and for the second sum, we use equation (A75) and obtain 

^i-m,k-m-n ~ ^N(i-m,k-m-n) ^ ( x i-k+n-l> y i-m-l) 


• (A96) 


where N(i,k) and Nj(ij) are defined in equation (A63) and the following equation. The 
variable subscripts were found from equation (A20) with i-m replacing i and k-m~n replacing 
k. 


For i > i c , we obtain from equation (A37) 


i-i c +k-2 


i-i c +k-3 i-i c +k-3-n 


C p " 2 C i c -m,l + 2 2 

m=0 c n=0 m=0 


^i c -m,i-i c +k-l -m-n 


i-ic- 1 


k-2 

H c -1 

I 

^i-m, 1 

+ Z 

2 „ C i 

m-0 


n=Q 

m=0 

i-i c -2 

i-i c -2- 

-n 


Z 

I 

C- 

n-m- 

-n— 1,2+m 

n=0 

m=0 




(A97) 
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In the first two sums, the C,- ^ terms are computed according to equations (A95) and (A96), 
respectively, for those panels lying in rows which terminate in a subsonic leading edge panel 
outside the zone of influence of the kth panel in the ith column. Separating these panels 
yields 


i i-i c +k-2 

k-2 i-V*- 3 ' 

2 

^i„-m,l + ^ ^ 

m=i-i c 

c n=0 m=0 

i-i c -i 

k-2 Hr 1 


+ X ^i-m,l + ^ s ^i-m,k+m-n 
m=0 n=0 m=0 


i'ic ^ 
+ 2 
m=0 


i-i c +k-3 



+ 2 

n=k-l. 


i-L+k-3-n 

2 c- • - , 

_ nu-m,i-u+k- 1 -m-n 
m-0 c t- 


i— i c — 2 i-i c -2-n 

S „ C i-m-n-l ,2+m (A98) 

n=0 m=0 


where the quantities outside the curly brackets require correction terms from subsonic 
leading edge panels. 

The pressure coefficient contribution for the first siim is computed from equation (A95) 
with i set equation equal to i c . For the second term we use equation (A75) and obtain 

^i c -m,i-i c +k = x N(i c -m,i-i c +k-l-m-n)^( x 2i c -i-k+n+2, y i c -m-l) (A99) 

where the variable subscripts were obtained from equation (A20) with i replaced by i c -m 
and k by i-i c +k-l-m-n. For the third term, we obtain from equations (A30) and (A81) 


^ 

i-m, 1 ^ X Nj(i-mj) S j ( x 2i c -i+m-l» y i-m-lj 


(A100) 


For i - i c > m > 0, we must add the correction given in equation (A85). Hence, we obtain 

VeJ_ 4 

C i c -m,l ” b X N i (i c -m,n) g n ( x i c -m-l , y i c -m-l) 


7rB [ XN l(i c “ m » 2 ) +Si X Ni(i c - m ,4) C ( X i c -m, y i c -m) 


(A101) 
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For m = 0, we also must include the correction in equation (A86). We obtain 


f t( X V y O[ X Nl(i c .l)' £ l S l X Nl(ic.2) +S l 


e l s l 2 x N,(i c ,4)J + B X Ni(i c ,n)8n (*i c -l, y i c -l) 

+ “il [ X Nl(i c >2) + s l X N!(i c ,4)j C (\, y i c ) 

for i - i c > m > 0, n > k-1, we add the correction given in equation (A85). This yields 
^i c ~m ,i-i c +k- 1 -m-n = 

nB X N(i c -m,i-i c +k-l-m-n)^ ( x 2i c -i-k+2+n , 
y i c -m-l)+ (sj/2) C (x 2ic _i_ k+2+n> yH c +k+2+n)J 


(A 102) 


(A103) 
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APPENDIX B 

INTEGRATION OF THE AERODYNAMIC INFLUENCE 
COEFFICIENTS FOR THE PLANAR MACH LINE PANELS 

B 1 BASIC INTEGRALS FOR SUPERSONIC LEADING EDGE PANELS 


For the basic integrals of the supersonic leading edge panels, we consider the following three 
integrals: designated by fj, f 2, 13: 


fl (x,y) = 


Vy 

f - dr> 

0 Vx/ej +(y-7? 2 ) 


\/y/( x fei + y) 

= f dt 

0 


= SIN' 




= fi (x,y) 

Since (x + e^i) /ej + (y - sj) = x/ei +y, we also have 

w 


Similarly, for the second integral we write 

f , , A ^ 2 d, 

f 2 (x,y)= j ■ — 

0 vx/ej+y-7? z 2 


x/ej +y 


= — s/yx/e'i + — j — f l ^ x ’ y ^ 


(Bl) 


(B 2 ) 


(B 3 ) 


and finally, for the third, we use 


Vy - 


f 3 (x,y) = / v 2 Vx/q + y-r? 2 dr? 

0 

yVyx/n . (x^+y) 

• = — 7” + v — z — 1 h ^y) 


(B 4 ) 
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The functions f 2 (x,y) and f,(x,y) have the same property as shown for fj(x,y) in equation 
(B2). We also require the-following limits on the' integration: 

■v/x/ej + y 

/ dr) r 1 dt • 

0 Vx/ej +y ~7? 2 ' Vi-t 2 2 


V x / e i + y 9 

r n z dn 

o Vx/ej +y-t? 2 


it (x/e 3 + y) y 
4 


V x / e i + y 2 

/ _2 r— j- / frfri +y) 

0 77 Vx/ej + y- 7 ? di? 16 


This is also interpreted as 

fl (x,y) = 7r/2 , x <0 

f2 (x,y) = ff (x/ej +y)/4 , x<0 

2 

f 3 (x,y) = 7r(x/ej + y ) fl 6 , x < 0 

B2 BASIC INTEGRALS FOR SUBSONIC LEADING EDGE PANELS 


For the subsonic leading edge, we consider the following three basic integrals: 






0 v 'x/e 2 ■ 


Vy 


f 5 (x,y) = / Vx/e 2 -y + T? 2 d7? = v ' yx ^ 2 + ^ f 4 (x,y) 


(B5) 


(B6) 


f 6 to) - /*" ,V x/e2 - y + „ 2 , f 5 0 ,y) 

o 4eo 


(B7) 


since (x/e 2 - s) - (y - s) = x/e 2 - y, we obtain for the limits Vy -s 2 t0 \Zy”. 


f 6 ( x -y) - % ( x ” e 2 s 2»y -s 2) 

Similar relations hold for f 4 and fj. 
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B3 AERODYNAMIC INFLUENCE COEFFICIENTS FOR SUPERSONIC LEADING 

EDGE DOUBLET PANELS 


Consider now the downwash at points for which- x > 0, y >-sj. In the panel the doublet 
distribution is of the form 

p(x,y) = (x+ejy) (a + bx + cy+ dxy) 
and the derivative with respect to y is 

/iy = (x+ejy) (c + dx) + ej (a + bx + cy + dxy) 

On the leading edge x = -ejy we obtain 

Py = ei £ a +(c-ei b )y-dei y 2 ] 


Also 


M xy = ( c + e l b) + 2d(x + e|y) 


Substituting equations (B8) and (B9) into equation (A5) and using the expression for the 
downwash. 


W 


= j[a gj (x,y) + b go (x,y) + c g 3 (x.y) + d g 4 (x,y)J 


leads to 


Mv/e7 

W = — - — - { 
v 


yr. 


t * r ti 

J I a + (c-e, b)(y- 77 2 )-de, (y-7r)~J dn/Vx/e, +(y-r? 2 ) 

[Vy-si 


yr 


v J r n 

2 J c + ej b + 4d (x + e] y-e]Tj 2 )/3 V x /e ] +(7- n 2 ) 

Vy^j L J • 


di7 


VT 


-2y/W\ J [ c + e, b+^p +2ej d(y-r? 2 )] 
v/y^sf 


dr] 


Applying the identities 


.r , — : — t ( x / e i +-y)-v 2 
Vx/ej + y - v 2 = 7 , 7 

Vx/ej + y-77 


t? 4 „ 2 / o . (x/ei +y)7? 2 

Vx/q+y-n 2 Vx/£ 1 + Y-* 


(B8) 

(B9) 

(BIO) 
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yields 


I yr ■ . . „ 

_ { J a+(c-6j b) y-dej y 2 dr?/ Vx/ej+y^" 

*r U/y - sj L . -* 


r vT 

- J (c - ej b - 2d ej y) t? 2 / v^x/ej +y - tj 2 
\/y-sj 

Vy“ r ' ' 

+ 2 X- [ C + € 1 b + 4d ( X + e , y) / 3 ] [( x /« 1 + V ) - V 2 ] *? / Vx/e ! + y - „ 2 


Vy-sj 


+ d , 


Vy~ 


vT 


/ V 2 \/x/e, + y -17 2 di? - (x/e, +y) / t? 2 djj/V x/e, +y-77 2 

Vy-sj 


Lvy^si 

2 \/x/e j £ c + ej b + 4d (x + ej y)/3 


(\/y-Vy-si) +4d e t S] V( x /e i)(y~si)/3 


V7 ; 

-8ejd J Tj^Vx/e] + y -17 2 d7? ► 
3 Vy-S! 


(Bll) 


'The result can be simplified even more by noting that, since x/ej + y = (x/e^ + sj) + (y-sj) , 

(x/ej+ y) J di? / Vx/ej +y-7j 2 =(x/e 1 +y) j^fj (x,y)-fj (x + ejSj, y -sj) J 

v/y-si 

= 2 [f 2 (x,y) - f 2 (x + y - sj)] 

+ v /yx/e 1 -V(y- Si ) (x/e 1 d-sf) 

Separating the gj functions according to-equation (BIO) yields 


gl (x,y) = f 1 (x,y) - fj (x + ejSj, y - sj) 


(B12) 
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g2(x.y) = e i (x » y)_f i ( x + e i s i> y " s i)J 

+ 3 [f 2 (x,y) - f 2 (x + ejSj, y - s^j 
+ 2y/y~- Si [Vx7ei - V x / e l +s l] 

g 3 (x,y) =y £fj (x,y)-fj (x + ejsj.y— S|)J 
+ f 2 (x,y)-f 2 (x + eis^y-sj) 

+ 2 N /y^ ^ V x / e l ~ V x / £ 1 +s l) 

§4 (x,y) = ej -y 2 [fi(x 3 y) -fj (x + ejs^y-s^J 

- (x/ej -y) |^f 2 (x J y)-f 2 (x + e 1 s 1 ,y-s 1 )j 
+ 8(x/ej+y) j^f 2 (x,y)-f 2 (x + e 1 s 1 ,y-s 1 ) 

+VV- S 1 {y*Je I - V x / e l +s l)J / 3 

- 5 ^f 3 (x,y) -f 3 (x + ejsj.y-s^j /3 
+ 4s iv '(y- s i)x/e 1 / 3 


(B13) 


(B14) 


(B15) 


Equations (B12), (B13), (B14), and (B15) hold when x > 0 and y > sj. Similar equations 
for x < 0, y < sj ; x < 0, y > s j ; and x >q, y < sj can be written down by setting the terms 
which become imaginary equal to zero and by assigning to fj(x,y).for x < 0 the limiting 
value as x goes to zero: i.e., 7r/2. Terms like 

V(x/ei)(y- Sl ) 

which are real for x < 0 and y < sj must be discarded also, because they represent solutions 
in the wrong nappe of the Mach cone. 
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B4 AERODYNAMIC INFLUENCE COEFFICIENTS FOR INTERIOR DOUBLET PANELS 
From equations (A8) and (B8) we obtain 


W 


= ^ / [a + 2b(x-{ 2 ) +2c(y-j, 2 ) 

+ 4d (x ~{ 2 )(y -T 7 2 )J d£d7? 


vy-*i V X T e I s l 


Integration yields 


w l2M 


|a + 4bx/3 + 4cy/3 + ^ dxyj( v /x’ — s/ x ~ e i s i) (V y -- W ~ s l) 


-2 sj (c + 4dx/3) N/y-s^V^ / 3 
-2ejSj (b + 4dy/3) V x - e l s l (>/y -\/y- s l ) / 3 

+ 4 dejS! 2 V( x - e I s l)(y- s l) / 9 

Defining the hj functions according to equation (A9) leads to 
hj (x,y) = 2 ^V x / e l “ V x / g l ” \/y - s lj 

h 2 (x,y) = 4xhj (x,y)/3 - 4e 1 s 1 V x / e l ~ s i (Vy - Vy~ s lj/ 3 
h 3 (x,y) = 4yhj (x,y)/3 - 4s j Vy-sj ^/x7q - V x / e l ~$i)/ 3 
h 4 (x,y) = (16/9) ^xyhj (x,y) - x Sj V(Y- s l) x / e l “ e l s l Y VY ( x / e l ~ s l) 

+ sj (x + ej y +ej sj/2) V(Y ~ s l)( x / e l ” s l) J 


(B 16 ) 


Equation (B16) holds for x > ejSj and y > Sj. The functions can easily be found for 
other ranges of x and y by taking the. real part, except when both y < s j and x < e j s j . The 
terms like 

VCy-Sj) (x/ej-sf) 

although real, must also be discarded since they are on the wrong side of the Mach cone. 
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B5 AERODYNAMIC INFLUENCE COEFFICIENTS FOR SUBSONIC 
LEADING EDGE DOUBLET PANELS 


Substituting equation (A 1.0) into.equation (A-l 2-)-yields 


W 


= | J £ a + (c + e 2 b) ( y - ?7 2 ) + de 2 (y - rj 2 ) 2 j d 77 / - y + V* 


VT /x-eiCy-T, 2 ) 

+ Ul^ 2 ) J J rc-e 2 b + 2d (x-e 2 y+ 

' Vy- s 2 V x - e 2 s 2 L 

This may be written, after one integration, 


-S 2 )] dfdTj 


W = 




yy 


- J fa + (c + e 2 b) ( x/e 2 ) + dx 2 /e 2 "| d7J ^ \A/e 2 -y + 

Vy^ 




+ J |"c + e 2 b + d (x + e 2 y)"| x^x/e 2 -y + t? 2 dr? 

Vy-sJ 



- d e 2 J t? 2 \/x/e 2 -y + t? 2 dr? 

Vy- S 2 


+ 2 


- e 2 b + 4d (x - e 2 y + e 2 t? 2 ) / 3 J \/ x / e 2 “ y + V* 


dr? 


Vy- s 2 

Vy 
/ 
vy- 


yy r /2x + e 2 s 2 A I 

_ 2 / c " e 2 b + 2d \ — 3 e 2 y + e 2 I V x /e 2 - s 2 dr? 

^7 c.2 L 


(B17) 


(B18) 
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Here we have made use of the identities 

y-i? 2 _ x / c 2 

\/x/e 2 - y + 17 2 v'x/e 2 -y + i?2 

(y-» 2 ) 2 ( x/£ 2) 2 . 2 s r — — , 

:: = -(x/e 2 + y-t? z )vx/e 2 -y + r? z 

y/x/e 2 -y + v 2 \/x/e 2 -y + 77 2 

Following the pattern of equation (A13) and equating coefficients of equation (B1 8) yields 
for the Sj functions 

Sj (x,y) = - f 4 (x,y) + f 4 (x - e 2 s 2 , y - s 2 ) 

5 2 (x,y) = x Sj (x,y) ~ e 2 J^f 5 (x,y) - f 5 (x - e 2 s 2> y - s 2 )] 

+ 2 e 2 V x / e 2 “ s 2 (VT- Vy ~ ^2 ) 

5 3 (x,y) = xS! (x,y)/e 2 + 3 [f 5 (x,y) - f 5 (x - e 2 s 2 , y - s 2 )J 

-2 Vx/e 2 - s 2 jyy"- Vy-S2 ) 

5 4 (x,y) = x 2 Sj (x,y)/e 2 

+(l lx - 5e 2 y) [f 5 (x,y) -f 5 (x- e 2 s 2 , y - s 2 )] 

+ (5g 2 /3) [f 6 (x,y) - f 6 (x - e 2 s 2 , y - s 2 )] 

+f- ( x - e 2 y) V x / e 2 ~ s 2 (vy-s2 ~ \/y~) 

4 e 2 s 2 i ; 

3 — Vy ( x / e 2 s 2) (B19) • 

Equation (B1 9) holds for x < e 2 s 2 and y > s 2 . For other ranges of the variables the relations 
are found by setting fy^fg equal to zero when one of the arguments is negative and 
dropping all radicals with negative arguments. The terms like 

V ( x / e 2 “ s 2) (y~ s 2) 

also must be discarded when x < e 2 s 2 and y < s 2 as discussed in the previous sections. 
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B6 AERODYNAMIC INFLUENCE COEFFICIENTS FOR SUPERSONIC 
LEADING EDGE SOURCE PANELS 

Substituting 

o(x,y) = a+bx+cy+dxy 


into equation (A69) yields 
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With the aid of the identities 


/ 


7]^ dr} 




+ y -v' 


-f 


(x/ej+ y)t ? 2 

-^•sA/ej + y-*r + / 

\ Vx/ej+y-r ? 2 


d 7 ? 


yV x / £ 1 + y-’J 2 drj -J 


- - f 3 (x,y) + (x/ej + y ) f 2 (x,y) 
(x/ej+y) dr? 


i? 2 di7 


-f - 

Vx/ej +y^7j 2 ,/x/ej + y^i? 2 


= f 2 (x,y) + v( x /q) y 

the preceding integrals can be expressed in terms of the basic integrals f j, f 2 , and f 3 
defined in equations (Bl), (B3), and (B4), respectively. Defining the functions gj, g 2 , g 3 , 
and §4 by means of equations (A70)', we obtain 

ii (x,y) = 2 (^p) [ f i - f i ( x + e i s i> y - s i)J 

g 2 (x,y) = €j y gj (x,y) + 2 (3 + ej) j^f 2 (x,y) - f 2 (x + ejSj, y - s^J 

g 3 (x,y) = y gj (x,y) - 2 jf 2 (x,y) - f 2 (x + ejsj, y - s x )J 

-2^/yxJel + 2V(y-si)(x/ej + s^j 

I 4 Cx,y) = - ej y 2 Ii (x,y) + 2 |- j^(3 - e t ) x/ej - (9 + ei ) yj Jf 2 (x,y) 

- f 2 ( X + ejsj, y - Sj)J+ (5 ej -3) Jf 3 (x,y) - f 3 (x + ejSj, y - Sj) J 
+ 2y (1 -ej) Vy x/ei 

+ 4 (x + y-sj) V (y -s l) x / e ! -2(2x- ei y + 3y-3s!) V(y-si)(x/e 1 +Si) 


/3 


(B20) 
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The preceding formulae hold for x > 0 and y > s j . Formulae for other regions of the flow 
may easily be written down by noting that • 


fj = 0 for y < 0. i-> -1 


and 


fj (x,y) = lim f j (x,y) = 7r/2 when x < 0 
x-> 0 

and by discarding all-radicals with negative radicands, including double radicals like 

V x y/<q 

when x and y both less than 0, even though they are real. 

B7 AERODYNAMIC INFLUENCE COEFFICIENTS FOR 
RECTANGULAR INTERIOR SOURCE PANELS 


For x > «2 s 2 anc * y > s 2> substituting 


a(x,y) = dxy 


into equation (A74) and using equation (A75) yields the following integrals to be evaluated 
for determining h(x,y) in equation (A75) 


y/y V3T 

h (x,y) = f f [x-S 2 + y-7j 2 ] 

Vy - s T v x_e i s i 


dfdi? 


Vy” V x_e i s i Vy~ s i' Vx" 

+ f. f ejSjdfdTjh f ' f S 1 dfdT? 

sjy - Sj 0 0 


Vx-ejsq 


This immediately integrates to 

h(x,y)=j [(x+y) y/vT-(x + y- s i) V x (y - s l) 


-(x + y-ejsj) Vy ( x " e l s T) 

+ (x + y-ejSj -Sj) V(y- S l) ( x “ e l s l)] 


(B21) 


The preceding equations hold for y > S 2 > x > ^ 2 - ^ ma y applied to other regions 
of the wing' by discarding all imaginary terms. As in the other influence coefficients terms 
like 

Vxy" 


when x < 0 and y < 0 are also discarded. 



B8 AERODYNAMIC INFLUENCE COEFFICIENTS FOR SUBSONIC 
LEADING EDGE SOURCE PANELS 


Substituting equation (A76) into equation (A80) yields the following relation for Cp 
the subsonic leading edge panel 


C p-S-iO-2) 


/ 


c(y-r? 2 ) +ej d(y-T? 2 ) 


Vy-S2 


2 N /x-e 2 (y-77 2 ) 


<377 


52 c s 2 + d e ? s 2 2 

0 


Vx-e 2 (y - r? 2 ) 


dr? 


Vy~ \/x-e 2 (y-f? 2 ) 

+ f jd (x - ^ 2 + y -7j 2 ) + c I dfdjj 

Vy~ s 2 V x -e 2 s 2 


Vy" V x - e 2 s 2 

+ (de 2 s 2 + c) f J dfdr? 

Vy-S2 O' 

Expressing Cp in the form of equation (A81) yields, for the functions S| and s 2 • 


Si&,y)=(-^) / Y (y-1 2 )^ tS2 / 

' 2 > ^x/« 2 -y + u z 2 y 


vy- s 2 


dr; 


s 2 

Vy - 

+ 2 / y/xfe 2 -y+v 2 d ?7 

Vy "- s 2 


>/x/e 2 -y + i? z 


S 2 fey) =(l-e 2 ) 


VT / 2 \2 o VT-S2 dTJ , 

f y-y*) <H? ^ 2 r dT? == 

x ^T 2 V /x / e 2-y + T ? 2 0 V x /e 2 -y + ^ 


Vy v^ x - e 2 (y-V 2 ) 

+7IJ / / (x-{ 2 + y-r? 2 )dfdj? 

Vy ~ s 2 V x “ e 2 s 2 
+ 2 e 2 s 2 V x /^2 - s 2 (vr- Vy - s 2 ) 


from 
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By using the relations 

n 

J-v 2 


x/e 2 * 


^/x/e 2 -y + Ti 2 x/x/e 2 - y + 17 2 


- y/xfe 2 - 


y+v" 


(y-v 2 Y 


(* 1 * 2 V 


\/x/e 2 -y + i? z \A/e 2 -y + 17^ 


- (x/e 2 + y - t? 2 ) \/x/e 2 - y + t ? 2 


the results of the integration may be expressed in terms of the functions f 4 , f 5 , and fg 
defined in equations (B 5 ), (B 6 ), and (B7), respectively. This leads to 


5 1 (x ’ y) = (^^).[( x/e 2) f 4. < x > y ) ~ ( x / e 2 - s 2-) f 4 C x " e 2 s 2’ y ” s 2>] 

+ ( el ^ [ f 5 ^ " f 5 ( x - e 2 s 2’ y - s 2) ] 

5 2 (x,y) = (1 - e 2 ) (x 2 /e 2 2 ) Jf 4 (x,y) - f 4 (x - e 2 s 2 , y - s 2 )J 

+ j(7 e 2 - 3) x/e 2 + (5e 2 + 3) y j (x,y) - f 5 (x - e 2 s 2 , y - s 2 ) j /3 

-(5e 2 + 3) j^fg (x,y) - fg (x - e 2 s 2s y-s 2 ) J/3 + ( 1 - e 2 ) s 2 2 f 4 (x - e 2 s 2 , y - 


-4(x + y-e 2 s 2 ) Vy(x/e 2 - s 2 ) /3 


+ 4(x + y-e 2 s 2 . + s 2 /2) v/(y~ s 2) ( x / e 2“ s 2)/3 - 

The preceding equations hold for y > s 2 and x > e 2 s 2 . Similar relations may be written 
down for other regions of the flow by following the rules in section B5 of appendix B. 


s 2> 

(B22) 
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B9 SUBSONIC LEADING EDGE CORRECTION TO THE AERODYNAMIC INFLUENCE 
COEFFICIENTS FOR SUPERSONIC LEADING EDGE SOURCE PANELS 
AND INTERIOR SOURCE PANELS 

Combining equations (A83) and (A84) and considering equation (A85), we see that the 
correction to g 2 , 14 , and h from the panel in those rows containing a subsonic leading edge 
in the zone of dependence of the point x,y is given by 


C(x,y)“ (1 -e 2 ) 


/ 


Vy - e 2 s 2 


(y-r? 2 ) dq 
Vx/e 2 - y + r? 


7 



-(2A/15-) 


W 

f f d$ dT ? 

\/y-s 2 >/x-e 2 (y- 7 ? 2 ) 


-(2V^) 



J d£dr? 

V x - e 2 s 2 


Performing the integration yields, after some simplification. 


C(x,y) = 


(l-e 2 ) (x/e 2 ) 


f 4 (x,y) -(x/e - s 2 ) f 4 ( x - e 2 s 2 



+ ( 1 + e 2> |^ f 5 < x >y) - % ( x - e 2 s 2> y “ s 2) J 

-2^(x/e 2 )y + 2 V(x/e 2 -s 2 ) (y-s 2 ) ( B23 > 


This relation holds for x>e 2 s 2 and y > s 2 . The expressions for other ranges of x and y are 
easily found by following the same procedure as for Sj and S 2 in section B5 of appendix B. 
The origin of the preceding equations corresponds to that for the subsonic leading edge 
panel in the same row as the panel for which the correction is to be applied. 
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Figure 2.— Mach Line Paneling on Planar Wing at M =\f2 
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(a) Control Point x,y Inside the Panel (b) Control Point x,y Outside the Panel 

and on y Characteristic Strip 


x,y Fixed point 

x Variables of integration 



x / 

(c) Control Pointx,y Outside the Panel (d) Control Point x,y Outside the Panel 

and on x Characteristic Strip and Between Characteristic Strips 


Figure 4.— Range of Integration for interior Mach Line Pane! and Analytic Continuation 
of Doublet Strength Outside of the Panel 
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Figure 5— Domain of Dependence of an Interior Pane! With Illustration of Upstream 
Characteristic Strips 
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Figure 6.— Control Point Locations and Continuity Conditions for Doublet Mach Line Paneling 
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Figure 7.— Analytic Continuation of the Doublet Distribution on the Wing Associated With 
a Mach Line Pane! Adjacent to Supersonic Leading Edge at M =\]~2 
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Figure 8.-Range of integration for Supersonic Leading Edge Mach Line Panel Analytic 
Continuation of Doublet Strength Outside of the Pane! 
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Figure 9.— Range of Integration for Subsonic Leading Edge Mach Line Panel and Analytic 
Continuation of Doublet Strength Outside of the Pane! 
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Figure 10.— Pressure Distributions on Portion of the Wing Behind Supersonic Leading Edge 
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Panel method 


-2,0 L, 


M = \/2 

t = maximum camber • 


Figure 1 1.— Pressure Distributions on Portion of the Wing Behind Subsonic Leading Edge 
Along the Line y = 1.0603 





Distance along line L in units of panel width 


Figure 1 3.— Discontinuity in Pressure Coefficient Along Panel Edges of Mach Line Doublet 
Panels in the Region Behind Subsonic Leading Edge 










Figure 1 5.— Pressure Distribution Along Line L Downstream of Special Mach Line With 
Fiat Plate Solution Superimposed on Doublet Panel Solution 
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Figure 16.— Pressure Distribution Along Line L Downstream of Special Mach Line With 
Flat Plate Solution Superimposed on Doublet Panel Solution 




Figure 17.— Pressure Distribution Along Line L Downstream of Special Mach Line With 
Flat Plate Solution Superimposed on Doublet Pane! Solution 




Figure 18.— Control Point Locations and Continuity Conditions for Source Mach Line 

Paneling 
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Figure 23 .— Region of Integration for Leading Edge Source Panel With the Zone of 
Dependence Terminated by a Subsonic Leading Edge 
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Figure 24 .— Region of Integration for Interior Source Pane! With the Zone of Dependence 
Terminated by a Subsonic Leading Edge 
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